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Abstract: We study in detail the general structure and further properties of the tree- level am- 
D ■ plitudes in the SU(N) nonlinear sigma model. We construct the flavor-ordered Feynman rules for 

!*i . various parameterizations of the SU(N) fields U(x), write down the Berends-Giele relations for 

the semi-on-shell currents and discuss their efficiency for the amplitude calculation in comparison 

£> \ with those of renormalizable theories. We also present an explicit form of the partial amplitudes 

up to ten external particles. It is well known that the standard BCFW recursive relations cannot 

^^ ' be used for reconstruction of the the on-shell amplitudes of effective theories like the SU(N) non- 
linear sigma model because of the inappropriate behavior of the deformed on-shell amplitudes at 

^T ■ infinity. We discuss possible generalization of the BCFW approach introducing "BCFW formula 

(f-} . with subtractions" and with help of Berends-Giele relations we prove particular scaling properties 
of the semi-on-shell amplitudes of the SU(N) nonlinear sigma model under specific shifts of the 
external momenta. These results allow us to define alternative deformation of the semi-on-shell 

S^ | amplitudes and derive BCFW-like recursion relations. These provide a systematic and effective 



tool for calculation of Goldstone bosons scattering amplitudes and it also shows the possible appli- 
cability of on-shell methods to effective field theories. We also use these BCFW-like relations for 
the investigation of the Adler zeroes and double soft limit of the semi-on-shell amplitudes. 
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1. Introduction 

The chiral nonlinear sigma model is a widely used tool for description of many phenomena in 
theoretical particle physics. It is based on a simple Lie Group G and the spontaneous symmetry 
breaking GxG — >■ G gives rise to massless excitations - Goldstone bosons. For instance, in the theory 
of strong interactions, the group G is SU(Nf) where Nf = 2,3 is a number of light quark flavors 
and Goldstone bosons are associated with the triplet of pions (for Nf = 2) or octet of pseudoscalar 
mesons it, K and r/ (for Nf = 3). The interactions of these degrees of freedom dominate the 
hadronic world at low energies. In this context, the leading order nonlinear U(3) x U(3) chiral 
invariant effective Lagrangian, the kinetic part of which corresponds to the chiral nonlinear U(3) 
sigma model, was constructed in the late sixties by Cronin [1] while the SU(2) case was studied by 
Weinberg [2,3], Brown [4] and Chang and Giirsey [5]. Further generalization lead to the invention 
of Chiral Perturbation Theory as a low energy effective theory of Quantum Chromodynamics by 
Weinberg [6] and by Gasser and Leutwyler [7], [8]. Chiral Perturbation Theory became a very 
useful tool for the investigation of the low energy hadron physics. 

The focus of this paper is on scattering amplitudes of Goldstone bosons within the SU(N) 
nonlinear sigma model described by the leading order Lagrangian. In principle, the standard Feyn- 
man diagram approach allows us to calculate arbitrary amplitude. Because the model is effective, 
and the Lagrangian contains an infinite tower of terms the calculation becomes very complicated 
for amplitudes of many external Goldstone bosons even at tree-level. It would be therefore desir- 
able to find alternative non-diagrammatic methods which could save the computational effort and 
provide us with a tool to get the amplitudes more efficiently. In the past an attempt to formulate 
the calculation of the tree-level without any reference to the Lagrangian was made by Susskind 
and Frye [9]. They postulated recursive procedure for pion amplitudes based on certain algebraic 
duality assumptions supplemented with the requirement of Adler zero condition which should have 
to be satisfied separately for group-factor free kinematical functions recently known as the partial 
or stripped amplitudes. Such a condition had been proven in the special case of pion amplitudes 
described by the SU{2) nonlinear sigma model by Osborn [10]. In [9] the authors successively 
calculated the amplitudes up to eight pions and showed that these results are equivalent to the 
diagrammatic calculation based on the SU(2) nonlinear sigma model. The full equivalence for all 
amplitudes has been proven by Ellis and Renner in [11]. 

Over the past two decades there has been a huge progress in understanding scattering am- 
plitudes using on-shell methods (for a review see e.g. [12-15]). They do not use explicitly the 
Lagrangian description of the theory and all on-shell quantities are calculated using on-shell data 
only with no access to off-shell physics (unlike virtual particles in Feynman diagrams). This has 
lead to many new theoretical tools (e.g. unitary methods [16,17], BCFW recursion relations for 
tree- level amplitudes [18, 19] and the loop integrand [20]) as well as practical applications of on- 
shell methods to LHC processes (for recent results of the next-to-leading order QCD corrections for 
W + 4-jets see [21]). Most of the recent theoretical developments have been driven by an intensive 
exploration of A^ = 4 super Yang-Mills theory in the planar limit both at weak and strong couplings 
(see e.g. [22-33]). 



There have been several attempts to extend some of these methods to other theories. The 
most natural starting point are the recursion relations for on-shell tree-level amplitudes, originally 
found by Britto, Cachazo, Feng and Witten for Yang-Mills theory [18], [19] and later also for 
gravity [34] , [35] . The main idea is to perform a complex shift on external momenta and reconstruct 
the amplitude recursively using analytic properties of the S-matrix. More recently, this recursive 
approach was extended to Yang-Mills and gravity theories coupled to matter, as well as more 
general class of renormalizable theories [36] . 

In this paper, we find the new recursion relations for all on-shell tree-level amplitudes of 
Goldstone bosons within SU(N) nonlinear sigma model. This shows that on-shell methods can be 
applied also for effective field theories and it gives new computational tool in this model. Using 
these recursion relations we are also able to prove more properties of tree-level amplitudes that are 
invisible in the Feynman diagram approach. 

The paper is organized as follows: In section 2 we discuss SU(N) nonlinear sigma model, 
introduce stripped amplitudes and using minimal parametrization (the convenient properties of 
which has been discussed in [11]) we calculate tree-level amplitudes up to 10 points. In section 3 we 
review BCFW recursion relations and their generalization to theories that do not vanish at infinity 
at large momentum shift. Section 4 is the main part of the paper, we first introduce semi-on-shell 
amplitudes, ie. amplitudes with n — 1 on-shell and one off-shell external legs. Then we prove scaling 
properties under particular momentum shifts which allows us to construct BCFW-like recursion 
relations. Finally, we show explicit 6pt example. In section 5 we use previous results to prove 
Adler zeroes and double-soft limit formula for stripped amplitudes. Additional results and technical 
details are postponed to appendices: In Appendix A, we describe the general parametrization of 
the SU(N) nonlinear sigma model. In Appendix B we give the results of the amplitudes up to 
lOp. Appendix C is devoted to the counting of flavor-ordered Feynman graphs needed for the 
calculations of the amplitudes in nonlinear sigma models and other theories. In Appendix D we 
present additional scaling properties of the semi-on-shell amplitudes. In Appendix E, we study the 
double soft-limit for more general class of spontaneously broken theories for complete (not stripped) 
amplitudes. 

2. Nonlinear sigma model 

2.1 Leading order Lagrangian 

Let us first assume a most general case of the principal chiral nonlinear sigma model based on a 
simple compact Lie group G. Such a model corresponds to the spontaneous symmetry breaking of 
the chiral group Gl x Gr where Gl,r = G to its diagonal subgroup Gy = G, i.e. to the subgroup 
of the elements h = (<7l, <7i?) where gi = gR- The vacuum little group Gy is invariant with respect 
to the involutive automorphism (<7l,<7.r) -^ (9r,9l) an d the homogeneous space Gl x G^/Gy is 
a symmetric space which is isomorphic to the group space G. A canonical realization of such an 
isomorphism is via restriction of the mapping 

(g L ,gR)^g R gl 1 = u (2.1) 



(which is constant on the right cosets of Gy in Gl x Gr) to Gl x Gr/Gv. Provided we induce the 
action of the chiral group on Gl X Gr/Gv by means of the left multiplication, the transformation 
of U under general element (Vl, Vr) of the chiral group is linear 

U^V R UV L \ (2.2) 

This can be used to construct the most general chiral invariant leading order effective Lagrangian 
in general number d of space-time dimensions describing the dynamics of the Goldstone bosons 
corresponding to the spontaneous symmetry breaking Gl x Gr — > Gy as 

c (2) = ^-(d^ud^u- 1 ) = —((U^dnUXU-^U)), (2.3) 

where F is a constant 1 with the canonical dimension d/2 — 1. Here and in what follows we use 
the notation (•) = Tr(-) and the trace is taken in the defining representation of G. The overall 
normalization factor is dictated by the form of the parametrization of the matrix U in terms of the 
Goldstone boson fields 4> a which we write for the purposes of this subsection 2 as 

U = exp (V2jA (2.4) 

where <f> = (p a t a and t a , a = 1, . . . , dim G are generators of G satisfying 

(t a t b ) = 5 ab (2.5) 

[t a ,t b ] = iV2f abc t c . (2.6) 

Here f abc are totally antisymmetric structure constants of the group G. According to (2.2), the 
fields (p a transform linearly under the little group Gy as the vector in the adjoint representation of 
G while the general chiral transformations of (p a are nonlinear. 

The Lagrangian CP 1 ' can be rewritten in terms of the Goldstone boson fields as follows. We 

, exp (- v / 2iAd(0)) - 1 1 exp (-|A/,) - 1 , N 

U » U Ad(0) W~ ^2 Dj, a4> [ZJ) 
where 

AAMdpt = [0, d^} = V2t a Dfd^ b ^V^t-D^- d<f>, (2.8) 

the matrix D^ b is given as 

Df = -i/ c< V (2.9) 

and the dot means contraction of the indices in the adjoint representation. Inserting this in (2.3) 
we get finally 

d , _ ^ ■ i^p ^ a, = -^ • (| fcff (I)" ^ . a,. ( 2 ,o, 



The decay constant of the Goldstone bosons. 
2 In what follows we will use also more general parametrization of U. 



2.2 General properties of the tree-level scattering amplitudes 

Note that, the only group factors which enter the interaction vertices are the structure constants 
f abc . In any tree Feynman diagram each f abc is contracted either with another structure constant 
within the same vertex or via propagator factor 5 with some structure constant entering next 
vertex. Therefore, using the standard argumentation for a general tree graph [12], i.e. expressing 
any / as a trace / = — (i[t a , t ]t c }/\/2 and then successively using the relations like f cde t c = 
— i[t d ,t e ]/ \pl in order to replace the contracted structure constants with the commutators of the 
generators inside the single trace, we can prove that any tree level on-shell amplitude has a simple 
group structure, namely 

M a ^- a -( Pl , P2 ,...,Pn)= Yl (t a ° (1) t a ° (2) ---t a ^ n, }M a ( Pl ,..., Pn ). (2.11) 

Here all the momenta treated as incoming and the sum is taken over the permutation of the n 
indices 1,2, ... ,n modulo cyclic permutations. As a consequence of the cyclicity of the trace we 
get 

M ff (pi,P2 ■ ■ ■ ,Pn) = M a (p 2 , . . . ,p n ,Pl) (2.12) 

Due to the Bose symmetry, the kinematical factors M. a {p\, , ■ ■ ■ ,p n ) has to satisfy 

M a0 p{pi, ,...,p n )= M<j{Pp(i),Pp(2),. • • ,P p (n)) (2-13) 

(where a o p is a composition of permutations) and therefore 

M ff (pi, ...,p n )= -M(p CT (l),p CT ( 2 ), • • • ,Pa(n)) (2-14) 

where we have denoted M = .Mid (here id is identical permutation). The amplitudes M{p\, . . . ,p n ) 
are called the stripped or partial amplitudes. Note that the same arguments can be used also for 
the Feynman rules for the interaction vertices, the general form of which can be written as 

V; ia2 - a "(pi, P 2, ■ ■ ■ ,Pn) = £ <tMW> . . . t a ^)V n (p a(1) ,p am ,- ■ ■ ,Pa(n))- (2-15) 

u£Sn/Z n 



After some algebra we get explicitly (see Appendix A for details) V2 n +i(pi, ■ ■ ■ ,P2n+i) = and 

V2n^...,p 2n ) = ^ (^) n_1 eW- 1 ( 2 ;: 1 2 ) !>■*+*)■ < 2 - 16 ) 

Let us note that besides (2.3), (2.4) we need not to use any algebraic relations specific for the 
concrete group G when deriving this formula and it is therefore valid for general G. In the general 
case we can therefore define the stripped amplitudes and stripped vertices, however, their relation 
is not straightforward and may depend on the group G. In what follows we will concentrate on the 
case G = SU(N). 



2.3 Tree-level amplitudes for G = SU(N) 

2.3.1 Flavor ordered Feynman rules 

The standard way of calculation of the tree-level amplitudes A^ ai '" a "(pi, . . . ,p n ) is to evaluate the 
contributions of all tree Feynman graphs with n external legs build form the complete vertices 
(2.15) and propagators A a b = i5 a b/p 2 . This includes rather tedious group algebra which is specific 
for each group G. In the special case of G = SU(N) the calculations can be further simplified. 
Because we have the completeness relations for the generators t a in the form 

N 2 -l 

£ (Xt a )(t a Y) = (XY) - jj(X)(Y), (2.17) 

a=l 

we can simply merge the traces from the vertices of any tree Feynman graphs in one single trace 
preserving at the same time the order of the generators t Qj inside the trace. Note that the "disco- 
nencted" 1/N terms have to cancel in the sum in order to produce the single trace in (2. II) 3 . This 
enables us to formulate simple "flavor ordered Feynman rules" directly for the stripped amplitudes 
M completely in terms of the stripped vertices V n . The general recipe is exactly the same as in the 
more familiar case of SU(N) Yang-Mills theory, i.e. the tree graphs built form the stripped vertices 
and propagators are decorated with cyclically ordered external momenta and the corresponding 
ordering of the momenta inside the stripped vertices are kept. 

Let us note that such a simple way of the calculation of the stripped amplitudes might not 
be possible for general group G. For instance for G = SO(N) we have the following completeness 
relations 

AT(JV-l)/2 

£ (Xn(t a Y) = -((XY)-(XY T )) (2.18) 

a=l 

the second term of which reverses the order of the generators in the merged vertex and the afore- 
mentioned simple argumentation leading to the flavor ordered Feynman rules has to be modified. 
The SU(N) case has also another useful feature. As a consequence of the completeness relations 
(2.17) for the group generators of SU(N) and the analogous relation 

N 2 -l 

£ (Xt a Yt a ) = (X)(Y) - -(XY) (2.19) 

a=l 

it can be proved [12] that the traces (t a °mt a °W . . . t a °^)) and (t a pd)£ a p(2) . . . £ a p(-)) are orthogonal 
in the leading order of N in the sense that 



Y (tMi)tMa) ...£<Mn))(t a P(D£ a P(2) ...t fl p(n))* =N n - 2 (N 2 - 1) (s ff p + o( — ] j (2. 



20) 



(21,(22, ...,an 



As we shall see in what follows, this fact can be understood as a consequence of the decoupling of the (7(1) 
Goldstone boson in the nonlinear U(N) sigma model. 



where b ap = 1 for p = a modulo cyclic permutation and zero otherwise. This relation is enough to 
uniquely determine the coefficients T a in the general expansion of the form 

j-aia 2 ...a n = V^ (£<Mi)£M2) . . . £<M«))7^, (2.21) 

aeSn/Z n 

(provided the coefficients T a are N— independent) as the leading in N terms of the "scalar product" 

^ r a ia2 ...a n ^a CT(1)t a CT(2) _ _ _ f a CT(n) y = jyn-2^2 _ y (j^ + q f J_ \ \ ^^ 

ai,a 2 ,...,a n ^ \ / / 

Because the stripped amplitudes and vertices by construction do not depend on N, the coefficients 
at the individual traces in the representation (2.11) are unique a therefore the stripped amplitudes 
and vertices are unique. 

2.3.2 Dependence on the parametrization 

Up to now we have identified the Goldstone boson fields (ft a using the exponential parametrization 
(2.4) of the group elements U(4> a ). However, according the equivalence theorem, the amplitudes 
J\/l a i a 2--- a n(p 1 ^p 2 ^ _ _ _ ^p n ^ are the same for any other parametrization U{<p a ) where 

4> a = <J) a + F a {(j)) (2.23) 

where F a (<p) = 0{(f) 2 ) is at least quadratic in the fields <fi. Therefore, according to the aforemen- 
tioned uniqueness, the stripped amplitudes for the nonlinear SU(N) sigma model do not depend 
on the parametrization. Note, however, that this is not true for the stripped vertices which do 
depend on the parametrization because the complete vertices V n x 2 " n (pi,P2, ■ ■ ■ ,Pn) do. 

As far as the on-shell tree- level amplitudes are concerned, in various calculations we are thus free 
to use the most suitable parametrization and consequently the most useful form of the corresponding 
stripped vertices for a given purpose. We shall often take advantage of this freedom in what follows. 

A wide class of parameterizations for the chiral nonlinear sigma model with G = U(N) and 
G = SU(N) has been discussed in [1]. The general form of such a parameterizations reads 

U = Y,ak[V2j(t>\ (2.24) 

where clq = a\ = 1 and the remaining real coefficients a^ are constrained by the requirement 
UU + = 1. The exponential parametrization (2.4) corresponds to the choice a n = \jn\. In fact, as 
was proved in [1], for SU(N) nonlinear sigma model with N > 2, the exponential parametrization 
is the only admissible choice within the above class of parameterizations (2.24) compatible with 
the nonlinearly realized symmetry with respect to the SU(N) chiral transformations (2.2). On 
the other hand, for SU(2) and for the extended chiral group G = U(N) with arbitrary N, the 
parameterizations of the form (2.24) represent an infinite-parametric class. The more detailed 
discussion can be found in Appendix A. 



2.3.3 Interrelation of the cases G = U(N) and G = SU(N) 

Let us note, that the SU{N) and U{N) chiral nonlinear sigma models are tightly related. Within 
the exponential parametrization we can write in the U(N) case 



[/ = exp \-\— <t>° \U (2.25) 




where U 6 SU(N) and <p° is the additional U{\) Goldstone boson corresponding to the U(l) 
generator t° = 1/yN. We get then 



U- l d^U = j\jjjd^ Q + U + dfi (2.26) 

and as a consequence, 

1 F 2 

£( 2 ) = -dcjP ■ d<j>° + —(d^Ud^U- 1 ). (2.27) 

Therefore cfP completely decouples. This means that for the on-shell amplitudes in this model 

M aia2 - an ( Pl ,p 2 ,...,Pn) = (2.28) 

whenever at least one a,j = 0. Note that this statement does not depend on the parametrization. 
We can therefore reproduce the on-shell amplitudes of the SU(N) chiral nonlinear sigma model 
from that of the U(N) one simply by assigning to the indices a« the values corresponding the SU(N) 
Goldstone bosons. Keeping this in mind, in what follows we will freely switch between the U(N) 
and SU(N) case and use the general parameterizations (2.24) also in the context of the SU(N) 
chiral nonlinear sigma model. 

The fact that the U(l) Goldstone boson decouples gives also a nice physical explanation why 
the "disconnected" 1/N term can be omitted in the relation (2.17) when summing over virtual 
states in the treedevel Feynman graphs for the SU(N) nonlinear sigma model. This term can be 
interpreted as the subtraction of the extra U(l) virtual state contained in the first "connected" 
part. However, because this state decouples, no such correction is in fact needed. 

The decoupling of the U(l) Goldstone boson is an effect analogous to the decoupling of the 
U(l) component of the gauge field in the case of the U(N) Yang-Mills theory. For the tree-level 
amplitudes (and the corresponding stripped amplitudes) we get as a consequence a set of identities 
constraining their form. For instance taking only one a,j = (say a\) in (2.28), we get the "dual 
Ward identity" (or the U(l) decoupling identity) 

M(pi,P2,P3,---,Pn) +M(p 2 ,Pl,P3,---,Pn) + ■■■ + M{p 2 ,P3, ■ ■ ■ ,Pl,Pn) = (2.29) 

exactly as in the Yang-Mills case (see e.g. [12] and references therein). 

2.4 Explicit examples of SU(N) on-shell amplitudes 

Using (2.11) we can reconstruct the complete amplitude .M ai '" 0n (pi, . . . ,p n ) just from a single 
stripped amplitude A4(pi, . . .p n ) which is given by the sum of Feynman diagrams with ordered 



external legs {1, 2, . . . n}. Though the aim of this paper is not to calculate scattering amplitudes 
using the Feynman diagram approach, in this section we provide explicit examples for diagrammatic 
calculation of the stripped 4pt and 6pt amplitudes of the chiral nonlinear SU(N) sigma model (the 
8pt and lOpt amplitudes we postpone to the Appendix B) as the reference result for the recursive 
formula given in section 4. 

We can easily see that the only poles in the stripped amplitude are of the form 1/sjj where 

j 
s i,j=P 2 i,j witn Pi,j = ^2Pk (2.30) 

k=i 

(Obviously Sjj = Sj+i^-i due to momentum conservation). The variables Sij are therefore well 
suited for presentation of the amplitudes. 

As we have discussed above, the SU(N) stripped amplitudes are essentially the same as those 
for the U(N) case and, as we have discussed above, they are independent on the parametrization 
of the unitary matrix U in (2.3). The most convenient one for diagrammatic calculation of on-shell 
scattering amplitudes is the minimal parametrization [11] 

U = V2^4> + \A-2^ = 1 + yfijt ~ 2 £ (^ ) k C„-i^ (2.31) 

fc=i ^ ' 

where C n are the Catalan numbers (A. 12). The stripped Feynman rules for vertices can be written 
in terms of Sij as follows (see Appendix A for details) 

ra-l 2n+2 



p2 J X Z-^t ^kCn-k-\ 2_j \i+2k+l 



V2n+2(Si,j) - ( ^2 ) o 7 .CfcCn-fc-1 / J Si,i+2fe+l (2.32) 



Note that within this parametrization the stripped vertices do not depend on the off-shellness of 
the momenta entering the vertex and when expressed in terms of the variables Sij they are identical 
taken both on-shell or off-shell. This rapidly speeds up the calculation, because there are no partial 
cancelations between the numerators and propagator denominators within the individual Feynman 
graphs and it allows us to find the final expressions for the amplitudes in very compact form. 

The four-point amplitude is directly given by the Feynman rule in the simple parametrization, 

2F 2 _M(l,2,3,4) = s li2 + S2>3 . (2.33) 

Note that for n-point amplitude X^fe=i Pfc = and this can be used to systematically eliminate p n 
or equivalently s. >n . 

The six-point amplitude is given by diagrams in Fig. 1. The explicit formula reads 

4F 4 A^(1,2,3,4,5,6) = 

(Sl,2 + S2,3)(S1,4 + S4,5) (S1,4 + S2,5)(S2,3 + S3,4) («1,2 + S2,5)(«3,4 + «4,5) 



«1,3 s 2,4 S 3j5 

+ (Sl,2 + Sl,4 + S 2 ,3 + «2,5 + «3,4 + «4,5) (2.34) 



10 




Figure 1: Graphical representation of the 6-point amplitude (2.34) with cycling tacitly assumed. 

This can be rewritten as 

4F*M(1, 2, 3, 4, 5, 6) = -\ ^ + ^^ + ^ + s lj2 + cycl , 
with 'cycl' defined for n-point amplitude as 



n-l 



A[Sij, . . . , s mj „J + cycl = 2_^ A[Si + k,j+k, ■ ■ ■ > s m+fc,n+fc] ) (2.35) 

fc=o 

which will quite considerably shorten the 8- and 10-point formulae. These are postponed to Ap- 
pendix B. 

3. Recursive methods for scattering amplitudes 

Feynman diagrams are completely universal way how to calculate scattering amplitudes in any 
theory (that has Lagrangian description). However, it is well-known that in many cases they are also 
very ineffective. Despite the expansion contains many diagrams each of them being a complicated 
function of external data, most terms vanish in the sum and the result is spectacularly simple. The 
most transparent example is Parke- Taylor formula [37] for all tree-level Maximal-Helicity- Violating 
amplitudes 4 . The simple structure of the result is totally invisible in the standard Feynman 
diagrams expansion. 

Several alternative approaches and methods have been discovered in last decades, let us mention 
e.g. the Berends-Giele recursive relations for the currents [38] and the more recent BCFW (Britto, 
Cachazo, Feng and Witten) recursion relations for on-shell tree-level amplitudes that reconstruct 
the result from its poles using simple Cauchy theorem [18], [19]. 

3.1 BCFW recursion relations 

For concreteness let us consider tree-level stripped on-shell amplitudes of n massless particles in 
SU(N) Yang-Mills theory ( "gluodynamics" ). 5 The partial amplitude Ai n is a gauge-invariant 
rational function of external momenta and additional quantum numbers h (helicities in case of 
gluons) 

M n = M n (pi,p 2 , ■ ■ -p n ; hi, h 2 , ■ ■ ■ h n ). (3.1) 



4 Scattering amplitudes of gluons where two of them have negative helicity and the other ones have positive helicity. 
5 The recursion relations can be also formulated for more general cases and also for massive particles. See [39] for 



more details. 
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The external momenta are generically complex but if we are interested in physical amplitudes we 
can set them to be real in the end. Let us pick two arbitrary indices i, j and perform following 
shift. 

Pi -+Pi(z) =pi + zq, pj ->• pj{z) = pj - zq (3.2) 

such that the momentum q is orthogonal to both pi and pj, ie. q 2 = (q • pi) = (q ■ pj) = and 
the shifted momenta remain on-shell. Let us note that such q can be found only for the case of 
spacetime dimensions d > 4. The amplitude becomes a meromorphic function A4 n (z) of complex 
parameter z with only simple poles. The original expression corresponds to z = 0. If M n (z) 
vanishes for z — > oo we can use the Cauchy theorem to reconstruct M n = .M n (0), 

1 f dz . , , . . , , , v^ Res (.M ra , zj.) 
0=—/ — Mn(z) = M n (0) + Y\ (3.3) 

2vri Jc(oo) z k Zk 

where C(oo) is closed contour at infinity. M n can be then expressed as 

M n = - £ ReS( - M "' Zfc) (3.4) 

k Zk 

where k is sum of all residues of Ai n (z) in the complex z-plane. Residues of M n {z) can be 
straightforwardly calculated for the following reason: the only poles of Ai n are p 2 b = where 
Pa,b = (Pa + Pa+i + • • -Pb)- The poles of Ai n {z) have still the same locations just shifted, namely 
Pa b( z ) = ^ where i G (a, a+1, ... 6) or j G (a, o+l, ... 6). If none of the indices i, j or both of them 
are in this range, the dependence on z in p a ,b{ z ) cancels and it is not pole in z anymore. It is easy 
to identify all locations of the corresponding poles z^. Suppose that particle i G (a, a+1, ... 6), 

o O s. n h 

Pa,b( z ) = {Pa + ■■■ Pi-l + {Pi + Zq) + p i+1 + . . . p b ) = =>• Z a)fe = - — ; (3.5) 

*{q ■ Pa,b) 

In the original amplitude M. n the residue on the pole p\ b = is given by unitarity: on the 

factorization channel with given helicity the amplitude factorizes into two sub-amplitudes, and 

therefore 

Res(.M n ,* a , b ) = y^M L (z a ^r h ^—^ K M h £ b (z a , b ) (3.6) 

f^ 2(q-p ab ) 

n ab 

where the summation over the helicities h ab of the one-particle intermediate state is taken. The 
"left" and "right" sub-amplitudes M^ h S b {z ab ) are 



l L,R 

! : ab ( „ . \ — \A. -l^ 



M L Sab (z a ^ b ) = M b - a +2(Pa, ■ ■ ■ ,Pi(z a , b ), . . .p b ,-p a ,b(Za,b);h a , . . . ,~h ab ) (3.7) 

M S R b (z a ,b) = M n -(p-a)(Pa,b(Za,b),Pb+l, ■ ■ ■ ,Pj(z a ,b), ■ ■ ■ ,Pa-l', h ab , ■ ■ ■ , h a -l). (3.8) 

The amplitude Ai n can be then written as 

M n = V Ml Kb (z a , b )^M h £ b (z a , b ) (3.9) 



,2 JVI R 

ah 



ab,h, ■- Pa < b 



It is convenient to choose i and j to be adjacent because it eliminates the number of factorization 
channels we have to consider. 
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3.2 Reconstruction formula with subtractions 

The BCFW recursion relations discussed above are very generic and applicable for a large class 
of theories. The main restriction is the requirement of large z behavior: Ai n (z) — > for z — > oo. 
However, this behavior is not guaranteed in general and there exist examples when it is broken no 
matter which pair of momenta pi and pj is chosen to be shifted. In such a case, an additional term 
(dubbed boundary term) is present on the right hand side of eq. (3.9). The boundary term, which is 
hard to obtain in general case, has been studied by various methods in the series of papers [40], [41] 
and [42], however no general solution is still available. Sometimes this problem can be cured by 
means of considering more general approach when all the external momenta p^ are deformed (such 
an all-line shift has been introduced in [43], see also [44]) 

Pk^ Pk(z) =Pk + zq k - (3.10) 

where z is a complex parameter and qt are appropriate vectors compatible with the requirements 
of the momentum conservation and on-shell constraint for Pk{z), ie. Pk' Qk = °\ = 0- The on-shell 
amplitude 

M n {z) = M n (pi(z),p 2 (z), • • • ,p n (z)) (3.11) 

become again meromorphic function of the variable z the only singularities of which are simple poles 
and the residue at these poles have the simple structure (3.6) dictated by unitarity. In some cases 
the desired behavior Ai n (z) — > for z — > oo can be achieved in this way. However, in general case 
the behavior of Ai n (z) for z — > oo is power-like with non-negative power of z. This fact requires 
some modification of the reconstruction procedure. 

This can be done as follows. Let us suppose that we have made any (linear) deformation of the 
external momenta p^ —> Pk{z) in such a way that the deformed amplitude M n {z) is a meromorphic 
function the only singularities of which are simple poles and let us assume the following asymptotic 
behavior 

M n {z)^z k (3.12) 

when z —> oo. Let us denote the poles of M n {z) as Zi, i = 1, 2, . . . n. Assume aj, j = 1, 2, . . . , k + 1 
to be complex numbers satisfying \a,j\ < R different form the poles z; L . Then we can write for z ^ aj 
inside the disc D(R) (i.e. inside the domain \z\ < R the boundary of which is a circle C{R) of the 
radius R) the following "k + 1 times subtracted Cauchy formula" (see Fig. 2) 



fc+i , 

f (../,>) w-z fJ^W-Uj 



27Ti rtm w — z - LJ - w — a 



fc+l fc+l .. , x fc+l n C(R) , .. s fc+1 

=A4n(*)iT— +r (aj0 tt -j_ + y Res (- M "-^ n^- (3-i3) 

i z - a j i a j - z , i ,, ■ a j - a i i Zi- z X1 ■ Zi- aj 

Here z\, z 2 , ■ ■ ■ , z nc(R) are the poles inside D(R) and Res (Ai n ; zi) are corresponding residues. In 
the limit R — >■ oo the integral vanishes due to (3.12) and D(oo) will contain all n poles. As a result 
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Figure 2: Illustration of the contour used for the derivation of the subtracted Cauchy formula (3.13) with 
k = 1 and nctm = 3. 

we get a reconstruction formula with k + 1 subtractions 



M n (z) 



y\ ReS (Mrcj Zj) i-r 2 

1=1 



fc+1 



fc+1 



n 

3=1 



+ ^M n (aj) J] 



Oj 



i=i 



l=i,¥j 



a, 



(3.14) 



This is the desired generalization of the usual prescription. In order to reconstruct the amplitude 
with the asymptotic behavior (3.12) from its pole structure, we need therefore along with the 
residues at the poles Zi (which are fixed by unitarity) also supplementary information, namely 
the k + 1 values A4 n (a,j) of the amplitude at the points a,j. Such a additional information is the 
weakest point of the relations (3.14): there exists no universal recipe how to get the values M n (a,j) 
for a general theory. This corresponds to the well known analogous situation of k + 1 subtracted 
dispersion relations, which allow to reconstruct a general amplitude from its discontinuities uniquely 
up to the fc + 1 generally unknown subtraction constants. Note that, provided we choose aj in such 
a way that M n (a,j) = (i.e. aj are the roots of the deformed amplitude Ai n (z)), we can reproduce 
the formula 

y^ Res {M n ; z t 



M n (z) 



i=\ 



Zi 



fc+1 

n 

3=1 



(3.15) 



first written in this context by Benincasa a Conde [45] and further discussed by Bo Feng, Yin Jia, 
Hui Luo a Mingxing Luo in [46]. 

4. BCFW-like relations for semi-on-shell amplitudes 

The straightforward application of the BCFW reconstruction procedure is not possible for the 
SU(N) nonlinear sigma model because the amplitudes M. n [z) do not have appropriate asymptotic 
behavior for z — > oo. The reason is that due to the derivative coupling of the Goldstone bosons 
the interaction vertices are quadratic in the momenta. Therefore after the BCFW shift the vertices 
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along the "hard" z— dependent line of the Feynman graph are in general linear in z and the linear 
large z behavior of the propagators cannot compensate for it. For instance, under the shift 6 (3.2) 
with i = 1, j = 2 we get for the 6pt amplitude (2.34) for z — > oo 

KA f \ o A9-P2,3)(si,4 + «4,5-«l,3) . fa ■ P2,5) {<l " P2,z) . {q " P2,d) («3,4 + «4,5 ~ ^3,5; 
yV46^J — — Z^ I 1 p r h 



«1,3 (g-P2,4) «3,5 

+0(A (4.1) 

and analogously A4 n (^) = O(^) f° r general 7 n. As discussed in the previous section, in order to 
reconstruct such an amplitude from its pole structure, it would be sufficient to know the values 
of M n {z) for two fixed values of z. However, such an information is difficult to gain solely from 
the Feynman graph analysis restricted only to the amplitudes A4 n . It is therefore useful to take 
into account also more flexible objects, namely the semi-on-shell amplitudes, which unlike the on- 
shell amplitudes depend on the parametrization of the matrix U and from which the on-shell 
amplitudes can be straightforwardly derived. As we would like to show in this section, appropriate 
choice of parametrization together with suitable way of BCFW-like deformation of the semi-on-shell 
amplitudes allows to substitute for the missing information on the amplitudes A4 n and to construct 
generalized BCFW-like relations for them. 

4.1 Semi-on-shell amplitudes and Berends-Giele relations 

The semi-on-shell amplitudes Jn ia2 '" an (pi,P2, ■ ■ ■ ,Pn) ( or currents in the terminology of the original 
paper [38], where they were introduced for QCD and more generally for the SU(N) Yang-Mills 
theory) can be defined in our case as the matrix elements of the Goldstone boson field (p a (0) 
between vacuum and the n Goldstone boson states \it ai {pi) ■ ■ ■ ir a ™(p n )) 

JZ aia2 - an (Pi,P2, ...,Pn) = W a (0)|7T ai (pi) • • • n a "(p n )}. (4.2) 

Here the momentum p n +i attached to <^> a (0) 

n 
Pn+1 = -^2PJ- (4.3) 

i=i 

is off-shell. Note that Jn ,aia2 '" an (pi,P2, ■ ■ ■ ,Pn) has a pole for p^+i = 0. 

In complete analogy with the on-shell amplitudes, at the tree level the right hand side of (4.2) 
can be expressed in terms of the flavor-stripped semi-on-shell amplitudes J n (pi,P2, ■ ■ ■ ,Pn) i n the 
form 

(0\r(0)\* ai (Pi)---K an (Pn))\trcc = Yl Tr(t a tMD ...^W)J„(p a(1) ,p a(2) ,..., Pa(n) ). (4.4) 



6 Under the all-line (anti)holomorphic BCFW shift the large z behavior is the same. Here we can use the general 
formulae derived in [44] which relate the number n of external particles, the sum H of their helicities and the overall 
dimension c of the couplings to the asymptotics of the amplitude under the all-line holomorphic (0(z a )) and anti- 
holomorphic (0(z s )) shift. These formulae reads 2s — 4 — n — c + H and 2a — 4 — n — c — H . In our case H — and 
the only coupling constant is -F -1 , therefore c = 2 — n, therefore in general case a — s — 1 independently on n. 

7 The general statement can be derived by induction from Brends-Giele recursive relations discussed in the next 
subsection. 
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J'l+1 



= E 




Figure 3: Graphical representation of the Berends-Giele recursive relations 

Let us note that, at higher orders in the loop expansion the group structure contains also multiple 
trace terms. We normalize the one particle states according to 

hip) = 1. (4.5) 

In this section the above semi-on-shell flavor-stripped amplitudes J n (pi,P2, ■ ■ ■ ,Pn) wm be the main 
subject of our interest. The on-shell stripped amplitudes M.{pi,P2, ■ ■ ■ ,Pn+i) can be extracted from 
them by means of the Lehmann-Symanzik-Zimmermann (LSZ) formulas 

M(pi,P2,...,Pn+l) = ~ lim pl +1 Jn(pi,P2,---,Pn)- (4.6) 

The main advantage of the semi-on-shell amplitudes J n (pi,P2, ■ ■ ■ ,Pn) (in what follows we also use 
short-hand notation J(l, 2, . . . , n)) is that they allow to abandon the Feynman diagram approach 
using appropriate recursive relation. The latter has been first formulated by Berends and Giele 
in the context of QCD [38] and proved to be very efficient for the calculation of the tree-level 
multi-gluon amplitudes. For the U(N) nonlinear sigma model the generalized recurrent relations 
of Berends-Giele type can be written in the form (see Fig. 3) 

n m 

J(l,2, ...,n) = -j- ^2 Yl iVm + 1 & 1 >h>Ph+hh>---P3rn-i+i,n,-Pi,n)Y[ J (Jk-i + h---Jk) ( 4 -7) 

P hn m=2 {jk} k =\ 

where the sum is over all splittings of the ordered set {1,2, ... ,n} into m non-empty ordered subsets 
{jk-i + 1, jk-i + 2, . . . ,j k }, (here j = and j m = n) 8 , V m+ i is the flavor-stripped Feynman rule 
for vertices with m + 1 external legs and p^k = Ylj=iPj as above. 

Let us note that, because the Lagrangian of the nonlinear sigma model includes infinite number 
of vertices with increasing number of fields, the above Berends-Giele relation for J n have to contain 



* Explicitly 

n — m + 1 n — m + 2 n — m+(m— 1) 

E- E E - E ■ 

{Jk} 3*1=1 j 2 =jl+l 3m-l=3m-2+l 
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n 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


t(2n+l) 


4 


12 


33 


88 


232 


609 


1 596 


4 180 


10 945 


28 656 


b(2n + l) 


5 


17 


50 


138 


370 


979 


2 575 


6 755 


17 700 


46 356 


/(2n + l) 


4 


21 


126 


818 


5 594 


39 693 


289 510 


2 157 150 


16 348 960 


125 642 146 


t 4 (2n+l) 


3 


6 


10 


15 


21 


28 


36 


45 


55 


66 


b 4 (2n+l) 


4 


10 


20 


35 


56 


84 


120 


165 


220 


286 


U(2n+l) 


3 


12 


55 


273 


1 428 


7 752 


43 263 


246 675 


1 430 715 


8 414 640 



Table 1: A comparison of the number t of the terms on the right hand side of the Bcrcnds-Giclc recursive 
relation with the total number b of terms needed for the Berends-Giele recursive calculation of the amplitude 
J(l, 2, . . . , 2n + 1) and with the total number / of flavor ordered Feynman graphs contributing to the same 
amplitude. In the last three row we compare these numbers with the analogous ones for the case of "0 4 
theory" . 

vertices up to n + 1 legs, i.e. much more terms than in the case of power-counting renormalizable 
theories like QCD where the number of vertices is finite 9 . This fact rather reduces the efficiency 
of these relation for the calculations of the amplitudes. We illustrate this in the Tab. 1, where 
the number of terms on the right hand side of the Berends-Giele relation (4.7) written for J2 n +i 
(denoted as t{2n + 1)) and the total number of terms necessary for the calculation of the same 
semi-on-shell amplitude using the Berends-Giele recursion (denoted as 6(2n + 1)) is compared with 
the total number f(2n + 1) of the flavor ordered Feynman graphs contributing to J211+1 an d with 
the same numbers valid for the theory with only quadrilinear vertices ("c/> 4 theory " ) denoted with 
subscript "4" . See Appendix C for more details and for derivation of the explicit formulae for 
these and other related cases. 

On the other hand, as we will see in what follows, the Berends-Giele relations can be used as 
a very suitable tool for the investigation of the general properties of the semi-on-shell amplitudes. 
Let us mention e.g. the following simple relations valid for J(l, 2, . . . , n) 



J(l,2,...,2n) = 
.7(1,2,... ,n) = J(n,n-l,...,2,l). 



(4.8) 
(4.9) 



These relation are valid independently on the field redefinition. However, as we shall see in what 
follows, some properties of the semi-on-shell amplitudes are not valid universally and are tightly 
related to a given parametrization. 

4.2 Cayley parametrization 

Unlike the on-shell amplitudes M ai " Mn (pi,P2, ■ ■ ■ ,Pn), which are physical observables and do not 
depend on the choice of the field variables provided the different choices are related by means 
of admissible (generally nonlinear) transformations, the concrete form of Jn' ai ' ari (pi,P2, ■ ■ ■ ,Pn) 
as well as the flavor-stripped amplitudes J n (pi,P2, ■ ■ ■ ,Pn) depends on the parametrization of the 



9 The number of terms on the right hand side of (4.7) grows exponentially with increasing n in contrast to the 
polynomial growths typical for the renormalizable theories. See Appendix C for details. 
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U(N) nonlinear sigma model. In what follows we will almost exclusively use the so called Cayley 
parameterizations 

where the Goldstone boson fields are arranged into the hermitian matrix (j) = <p a t a with t a being 
the U(N) generators. As described in Appendix A, representation (4.10) is a special member 
of a wide class of parameterizations suited for the construction of the flavor-stripped Feynman 
rules. The interrelation between the field </> and analogous field (j) of the more usual exponential 
parametrization U = exp f j?4>) is through the following admissible nonlinear field redefinition 

«/, = 2FtanQ^) = ^+C>(>). (4.11) 

As is shown in Appendix A, the flavor-stripped Feynman rules for vertices read in the Cayley 
parametrization 

vw = o 

r_n" / 1 \ 2n n 2n+2 f— i \ n f 1 \ 2n / n \ 2 

V2 n+ 2 = - [ - ¥ ^ r U) EE(w^+i) = i 2r (f) U> +1 ' (412) 

V y jr=0 i=l V y \i=0 / 

where we have used the momentum conservation in the last row. For the first non-trivial vertex V4 
we get 

V 4 = -^L( Pl+P3 f = __L( P2 +p4 )2 (4.13) 

and the first two non-trivial semi-on-shell amplitudes read in the Cayley parametrization 

J(1,2,3) = 2i^y (pi + prf ( 4 - 14 ) 



.7(1,2,3,4,5) 



1 



AF 4 p 2 



(Pi + P2 + P3 + P5)(Pl +P3F , {P1+P3+P4+PS)(P3+P5) 



+ 



(P1+P2+P3) 2 ' (P3 + P4 + P5) 2 

(4.15) 



(P1+P5) 2 (P2+P4) 2 f , , v 2 

H 7 i ! Y2 (P1+P3+P5) 



(P2 +P3+ PA) 



Let us illustrate explicitly the dependence of the semi-on-shell amplitudes on the parametrization. 
Using the exponential one we obtain different amplitude J(l, 2,3), namely 

j,, « on 1 (gl + P2? + (P2 + P3? ~ 2(gi + p 3 ) 2 ,, 1R v 

J(l,2,3) eX p=-g^2 ^2 • ( 4 - 16 ) 

However, both J(l,2,3) and J(l,2,3) exp give the same on-shell amplitude (2.33). 

In the next subsection we will prove additional useful properties of the semi-on-shell amplitudes. 
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4.3 Scaling properties of semi-on-shell amplitudes 

The Cayley parametrization is specific in the sense that the semi-on-shell amplitudes J n (pi, ■ ■ ■ ,p n ) 
in this parametrization obey simple scaling properties when some subset of the momenta pi are 
scaled pi — > tpi and the scaling parameter t is then send to zero. Here we will study two important 
scaling limits, corresponding to the case when all odd or all even on-shell momenta are scaled. As we 
shall see in the following section, these two scaling limits are the key ingredients for the construction 
of the BCFW-like relations for semi-on-shell amplitudes in the Cayley parametrization. 
We will prove that for n > 1 and t — >■ 

J2n+l(tpi,P2,tp 3 ,P4, ■ ■ ■ ,P2r,tp2r+l,P2r+2, ■ ■ ■ ,P2n,tp 2n +l) = 0(t 2 ) (4.17) 

and 

lim J 2n +l(Pl, tp2,P3, tp4, ■■■, tp2r,P2r+l,tp2r+2, • • • , tp 2 n,P2n+l) = 7^p2W ( 4 ' 18 ) 

The general proof of (4.17) and (4.18) is by induction. Let us first verify the base cases. While 
the second statement holds already for n = 1 

J3{Pl,tp2,P3) = T^7 — : T2^77^2> ( 4 - 19 

F \P\ +tp2+Psr ^ F 

the first one is not valid unless n = 2. Indeed 

Mtpi,P2,t P3 ) = -^- ^i-P3) - = 0(t). (4.20) 

2F 2 (pi • p 2 ) + (p 2 • P3) + t{Pl ■ P3) 

On the other hand, using the explicit form of J5 (cf. (4.15)) we get 

Js{tpi,P2,tpa,PA,tp5) = 0{t 2 ); (4.21) 

we can therefore proceed by induction starting at n = 2. 

Let us first prove the scaling property (4.17). Suppose, that (4.17, 4.18) holds for all n, where 
1 < n < n and write for the left hand side of (4.17) the Berends-Giele relation (4.7) expressing 
J 2 n+i i n terms of J2n+i with n < n. After the scaling p 2 k+i -^ tP2k+i, the t — > behavior of p 2n+2 
and Kn+i is O(t ) and 0(t r ) where r > respectively. The scaling of the remaining semi-on-shell 
amplitudes on the right hand side of (4.7) can be deduced from the induction hypothesis. Note 
that it depends on the number of the external on-shell legs of J{ji-i + 1, . . . ,ji) as well as on the 
parity of jj_i + 1, because the semi-on-shell amplitude with scaled even or odd momenta scales 
differently. Namely, according to the induction hypothesis, the scaling of these building blocks of 
the right hand side of (4.7) is as follows (see Fig. 4) 

J(j) = 1 = O(t ), J(2j-l,2j,2j + l) = 0(t), J(2j,...,2k) = O(t ), 

J(2j + 1, . . . , 2k + 1) = 0(t 2 ) for k - j > 1. (4.22) 

This implies, that those terms of Berends-Giele relations which are depicted in Fig. 5, i.e. those 
which contain at least one block J(2j + 1, . . . ,2/c + l) = 0(t 2 ) with k — j > 1 or at least two building 
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2j 




o(t°) 



2J-1 




2 J = 0(t) 



2j + l 




0(t° 



0(t 2 



2k + 1 



Figure 4: Scaling of the building blocks on the right hand hand of the Berends-Giele recursion relation 
according to the induction hypothesis when the odd momenta are scaled. 



2ji - 1 




2ji - 1 



0(t< 



2j m _i + 1 




0(t 2 



2n + l 



Figure 5: The terms on the right hand hand of the Berends-Giele recursion relation which are automatically 
0(t-) using the induction hypothesis when the odd momenta are scaled. 

blocs J(2j — 1, 2j, 2j + 1) are automatically 0(t 2 ). Therefore, the only dangerous terms on the right 
hand side of (4.7) are those without the buildings block of the type J(2j + 1, . . . , 2k + 1) = 0(t 2 ) 
with k — j > 1 and at the same time without (case I) or with just one (case II) building block 
J(2j — 1, 2j, 2j + 1) = 0(t) (see Fig. 6). To this terms the induction hypothesis cannot be applied 
directly. 

In the case I, the odd lines of the corresponding vertex V2m+2 are attached to J{2jk + 1) = 1 
and such a vertex is then proportional to the squared sum of the odd momenta tp2j k +i, (cf. (4.12)) 

V2m+2{tpi,P2,2j 1 ,tp 2 j 1 +l, ■ ■ ■ ,tp2n+l) ~ (*Pl + tP2ji+l H \~ tp2n+l) 2 (4.23) 

which means that it scales as 0(t 2 ). This is in fact the scaling of the complete contribution of the 
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case I 



case II 




2n + l 



2n + l 



V 2m +2 ~ (t'Pl + t P2n + l + ...+ t P2n+1 ) 2 = 0(t 2 ), V 2m +2 ~ (iP2.,-l + P2, + tp 2j +l + Efc iP2.,,+l) 2 = 0(t l ) 

Figure 6: Typical terms on the right hand hand of the Berends-Giele recursion relation to which the 
induction hypothesis (4.17) cannot be applied directly. In both cases, to all (case I) or to all but one (case 
II) odd lines of the vertex the blocks J\ are attached. In the case II, one building block J 3 is attached to 
remaining odd line. 

terms in the case I, because all the remaining building blocs are of the order 0(t°) for t — > 0. 

In the case II with exactly one building block Js(tp2j-i,P2jitp2j+i) = 0(t) (note that, it has 
to be attached to the odd line of the vertex V2m+2), all the other odd lines of V^m-^-i are attached to 
J(2jk + 1) = 1 and such a vertex is then proportional to the squared sum of the momenta tp2j k +i 
and the momentum of the line which is attached to Js(tp2j-i, P2j ,tp2j+i) , namely 



V> 



2m+2 



( tP2j-l + P2j + tP2j+l + Yl tPlJk+l j 



0(t). 



(4.24) 



Therefore the complete contribution of the dangerous terms in the case II is in fact 0(t 2 ) for t — > 
because both V2 m +2 and Jz{tP2j-i-,P2j-,'tP2j+i) scale as 0(t) and again all the remaining building 
blocks are of the order 0(t°) for t — >• 0. All the other "non-dangerous" terms on the right hand 
side of the Berends-Giele relations scale at least as 0(t 2 ), which finishes the proof of (4.17). 

Let us now prove (4.18), i.e. the case when all even momenta are scaled. Suppose validity of 
this relation for n < n and again write the Berends-Giele relation for the left hand side of (4.18). 
Thanks to the just proven statement (4.17), the terms on the right hand side of (4.7) with at least 
one building block J{j^ + 1, . . . ,jk+i) with odd jf. and jk+i — jk > 1 do not contribute in the limit 
t — > 0. Such a block can be attached only to the even line of the vertex V m +i. Therefore, the only 
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terms which can contribute in the limit t — > have the form depicted in Fig. 7, i.e. those with the 
building blocks Ji attached to all even lines of the vertex. 



2ji - 1 




2ji + 1 



2j2 - 1 



27„ + 1 



Figure 7: Typical terms on the right hand hand of the Berends-Giele recursion relation which contribute 
to (4.18). Here to all even lines of the vertex the blocks J\ are attached. 

According to the induction hypothesis and using the explicit form of V2k+2 this gives for t — > 

\k k + 1 1 I -i\k 



-iy 



11 (OF2Y7I-.7J-1-1 



(-1)' 



2kp2k 11 {2F^\n~ii-i-"^ 2 n F 2n 

where we denote jo = and j^+i = n + 1. Sum of all such contributions is 



n ( -nfc-i 

/ j / ^ 2 n F^ n 2 n F^ n / J \ k 

k=ll<ji<J2<—,Jh<n ' fc=l V 



( - ir '- ' 5 i:[:](-r 



2"n.p2r, 



(4.25) 



(4.26) 



which finishes the proof. 

Another independent scaling properties of the semi-on-shell amplitudes J2n+i can be proven 
using the same strategy. For instance, when all odd momenta and one additional even momentum 
(say p2r) are scaled, we get 

lim J 2 n+l(tpi,P2, tp3,P4, ■■■, tp2r-l,tp2r,tp 2 r+l, ■ ■ ■ ,P2N,tp2n+l) = (4.27) 

for n > 1. We postpone the proof to the Appendix D. 

Let us note that due to the homogeneity of J(l, 2, . . . , 2n + 1) we can rewrite the relations 
(4.17) and (4.18) as a statement on the asymptotic behavior of the scaled amplitudes for t — > oo, 
namely 

lim J 2n +l(tpi,P2, ■ ■ ■ ,P2n,tp2n+l) = hm J2n+l(j>l,t~ 1 p2, ■ ■ ■ ,t~ 1 p2n,P2n+l) = , 2/? 2\n ( 4>28 ) 



t— >co 



t— >oo 



and 



J2n+l(Pl,tp2, ■ ■ ■ ,tp2n,P2n+l) = ^2n+l(* Pl,P2, ■ ■ ■ ,P2n, t l P2n+l) = 0{f 



(4.29) 
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4.4 BCFW reconstruction 

As we have mentioned in the previous subsection, the standard BCFW-like deformation of the 
external momenta pi yields deformed amplitudes which behave as a non-negative power of z for 
z — >• oo. As a result, for the reconstruction of the amplitude from its pole structure we need to use 
the general reconstruction formula (3.14) for which additional information on the on-shell amplitude 
(its values at several points) is necessary. However, such an information is not at our disposal. We 
solve this problems by the following trick: we relax some demands placed on the usual BCFW- 
like deformation and allow more general ones for which either the reconstruction formula without 
subtractions can be applied or additional information on the deformed amplitudes is accessible. 
The momentum conservation cannot be evidently avoided, what remains is the on-shell condition 
of all the external momenta. It seems therefore to be natural to relax this constraint and instead of 
the on-shell amplitudes M.2n+2 to use the semi-on-shell amplitudes J2n+i, or the cut semi-on-shell 
amplitudes M.2 n +i defined as 

M 2n +1 (Pi, • • ■ ,P2n+l) = Pl,2n+1-Wl (PU ■ ■ ■ ,P2n+l) • (4.30) 

Motivated by the results of the previous section let us assume the following deformation of the 
semi-on-shell amplitude M2 n +i in the Cayley parametrization 

M 2n +l{z) = M 2n +l(Pl,ZP2,P3,ZP4, ■ ■ .,Zp 2r ,P2r+l,Zp2r+2, ■ ■ ■ ,ZP2n,P2n+l) (4.31) 

i.e. all even momenta are scaled by the complex parameter z and the odd momenta are not deformed 

P2k{z) = Zp2k, P2k+l(z) =P2k+l (4.32) 

Note that in contrast to the standard BCFW shift this deformation is possible for general number 
of space-time dimensions d. The physical amplitude corresponds to z = 1. For n = 1 we get 
explicitly 

M 3 (z) = j^( Pl - P3 ) (4.33) 

For general n let us denote the sums of all odd (even) momenta as 

n n 

P- = ^P2k+i, P+ = ^P2k ■ (4.34) 

k=0 fc=l 

Then in general case the function M2 n +i(z) has the following important properties: 

1. With generic fixed pi it is a meromorphic function of z with simple poles. 

2. The asymptotics of M2 n +\{z) can be deduced form the known properties of J2n+i, namely 
for n > 1 we get as a consequence of (4.29) 

M 2n +l{z) = {p+Z + p-) 2 J 2 n+l(pi, zp 2 , ■ ■ ■ , Zp 2n ,P2n+l) = 0(z°). (4.35) 
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3. For n > 1 we have according to known scaling property (4.18) of J2n+i 

limAf 2n+1 (z) = ^l^p 2 _ (4.36) 

The first two properties allows us to write for M2 n +i{z) the reconstruction formula with one 
subtraction, i.e. the relation (3.14) with k = 0. The third property is the key one for the complete 
reconstruction and determines both the "subtraction point " a± = and the "subtraction constant " 
M 2n+ i(ai) = p 2 _/(2F 2 ) n . The resulting formula reads 10 

*wi(*) = ^p- + E — jr^ — - ( 4 - 37 ) 

where the sum is over the poles zp of M2 n +i(z). The position of the poles is known and the 
corresponding residues can be determined recursively as in usual BCFW relations, however, there 
are some subtleties. 

The poles zp of M2 n +i(z) correspond to the vanishing denominators of the deformed propaga- 
tors pp{z) = 0, where 

p 2 P (z) = PiJ (z) 2 = 0, for 2<j-i < 2n (4.38) 

and where j — i is even; in this formula Pij(z) = zpf- + p~- with 

p tj = E P2k > Pij = E P2*+l' ( 4 ' 39 ) 

i<2k<j i<2k+l<j 

i.e. Pj ■ is a sum of all even (odd) momenta from the ordered set Pi,Pi+i, ■ ■ ■ ,pj~\,pj. Explicitly 
for j — i > 2 

1/2 



z?j = ^ ^— (4.40) 



,2l2 („ . ln2 



where G(a,b) = a b — (a ■ b) is the Gram determinant, which is nonzero for generic momenta 
Pi, . . . ,Pj. Therefore in the generic case f or j — i > 2 we deal with doublets of single poles. 

The case of three-particle poles corresponding to j — i = 2 has to be treated separately. In 
this case either pf- = or p~ ■ = (this sets in for pf- = p i+ i or for p~ ■ = p i+ i respectively; let us 
remind that pk are on-shell). In the first case we have only one pole 

_ (P2J-1 ■ P2j+l) UAU 

Z2j-l,2j+l — / ; r I 4 - 41 ) 

P2j-{P2j-l+P2j+l) 



1 Let us note, that we could write analogous reconstruction formula directly for the currents J2n+i as we did 
in [49]. In such a case we do not need any subtraction. The price to pay is that we get two more poles, the residues 
of which cannot be determined recursively from unitarity. Fortunately, the relation (4.29) and the residue theorem 
can be used in order to obtain the unknown residues in terms of the remaining ones. The resulting formula is fully 
equivalent to (4.37), however it is a little bit less elegant. 
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ptA*&) 




In A- 1 



j + 1 



Figure 8: Graphical representation of the right hand side of the relation (4.45). 



while in the second case we have apparently two poles 



'2j,2j+2 



Z 2j,2j+2 = Z 2j,2j+2 



P2J+1 ■ {P2j + P2J+2) 
(P2j ■ P2J+2) 



(4.42) 
(4.43) 



However z^ 2 ■ , 2 = cannot be a pole according to (4.36) and the corresponding residue has to be 



zero. 



The residues of the function M 2n+ i(2:) are dictated by unitarity and at the poles they factorize 
(see Fig. 8). Writing for j — i > 2 



(zpfj + P hj f = pf](z ~ «£•)(* ~ Z i,j) 



we get for j — i > 2 



Res [M 2n+ i,z i 



M^Hz^M^Hz*) 



(4.44) 



(4.45) 



-M7 " v +2 (z + - z~ ) 

Pi,j\ z i,j z i,j> 

where we denoted 

M L J \4j) = Mh+l-(j-i)(Pl(4)' • • • >Pi-l( Z i 3 )>PiA Z i 3 )>Pj+l( Z tj)i ■ ■ ■ ,P2n+l(zf; j )) (4.46) 

M^\zf.) = M^ i+1 { Pi {zt tj ), Vi+l {zt tj ), . . . >Pi (z£))- (4.47) 

Note that, while the amplitude M^' remains semi-on-shell, the amplitude M^' is fully on-shell, 
because the deformed momentum Pij(z) is on-shell for z = z i ■. 

The formula (4.45) is valid also for the three-particle pole Z2j,2j+2 given by (4.43). However 
the pole Z2.7-i,2j+l deserves a special remark because the corresponding residue is determined by 
the formula different from (4.45), namely 



Res(M 2n+1 ,z 2i _ li2j+1 ) = — L - \ \ h2j+l) g 

2 P2j-l,2j+l ■P2j-l,2j+l 



(Z2j-l,2j+l) 



(4.48) 
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where M L 3 R ' 3 (z2j~i,2j+i) are given by (4.46) and (4.47) with z i • replaced by Z2j-i : 2j+i- 

To summarize, we have found a closed system of recursive BCFW-like relations for the tree 
cut semi-on-shell amplitudes M2 n +i, which consists of the reconstruction formula (4.37), the pole 
positions (4.40), (4.41) and (4.43) and the residue formulae (4.45) and (4.48). Note that the initial 
condition for the recursion (4.33) can be understood as the special case of (4.37) for n = 1 because 
then there is no pole Zij with 2 < j — i < 2 and the sum of the residue contributions is empty. The 
physical amplitude M2 n +i(pi, ■ ■ ■ ,P2n+i) corresponds to z = 1 

M 2n+1 ( Pl , . . . , P2n+ i) = J^P- + L 73^ (449) 

As a final result we get then using (4.45), (4.48), (4.41), (4.43) and (4.44) 

M 2n+ l(pi,...,P2n + l) = 1 ^^P 2 -+Y. M L ) ^P)^T M R ) ^P)- ( 45 °) 

\ zr ) p Pp 

Note that there is an extra function Rp in contrast to the standard BCFW formula (3.9), namely 



Rp 



for Zp = Z 2 j,2j+2 



-1 



for z P = z 2j -i,2j+i (4.5i) 



1 1 Z S- for z P = zf 



*lj *J ',J 

For further convenience, we rewrite (4.50) with help of (4.33) in the following more explicit form 
M 2n +i {pi,... ,P2n+l) = (2F 2 ) n P - + 

4. V^ A/f (2i,2j+2) . v 1 P2j • P2J+2 
+ Z^ M i (22j,2j+2j-2 



i=i 
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'2 p2 

P2j,2j+2 r 

(2j-l,2j+l), x 1 P2j-l,2j+l'P2j-l,2j+l 



j;«f- 1A+1 '(z, 



- Z^ iW L ^2-,-l,2j+i;-2 ^2" 

j=l P2j-l,2j+l ^ 

2<i-i<2n Z i,j Z i,j V ^J \j P i,j Z i,j / 

(4.52) 

The on-shell amplitude is then 

M 2n (l,2,...,2n-l;2n) = - lim M 2n -i(l). (4.53) 

4.5 Explicit example of application of BCFW relations: 6pt amplitude 

As an illustration let us apply the BCFW-like recursive relations (4.37) to the amplitude M-(z) = 
Ms(pi, zp2,P3, zp4,ps). In this case we have three poles, all of them being three-particle, namely 

1 Sl ' 3 ft S2 ' 4 ^ 1 S3 ' 5 (A KA\ 

21,3 = 1 7 , 32,4 = U 7 , Z3,5 = l 7 (4-54) 

•51,2 + 52,3 V 52,3 + 53,4/ «3,4 + «4,5 
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where the variables Sij are given by (2. 30). The residues are given by the relations (4.45) for Z2,4 
and (4.48 ) for z\^ and 2:3,5. After simple algebra using the explicit form of the poles (4.54) we get 



Res(M 5 ,zi, 3 ) 

Zl,3 

Res (M 5 , ,23,5) _ 

2 3 ,5 ~ 4F 4 

Res(M 5 ,z 2 ,4) 1 



'1 - Zl,3) (S2,5 - «2,4 + S3,4 ~ S 3)5 ) - — =j (si j5 - Sl ;4 - S 4j5 ) 



4i?4 

1 



4i?4 



'1 - Z3,5)(S1,4 - Sl,3 + s 2,3 - S2,4) 



(Sl,5 - Si 4 + S2A ~ S 2 ,5) ■ 



4i?4 



(Sl,5 - «1,2 - S2,5) 



z 2A 4F 4 

Note that the potential unphysical poles Zij(pk) = have canceled completely. We have also 

,-1 Sl,2 + S2,3 /-, x-l S 3; 4 + S 4; 5 



(4.55) 



(1 " 21,3)" 



Sl,3 



-, a-^3,5)- 1 



S3,5 



(1 - -22,4) = 1 ! (4.56) 



«2,4 



These factors are responsible for setting of the physical poles in the resulting amplitude. After 
inserting this to the formula (4.49) we get for the individual contributions to the semi-on-shell 
amplitude in the Cayley parametrization 



Res(M 5 ,zi, 3 ) __ _L_ 

21,3(1 - 21,3) 4F 2 

Res(M 5 ,z 3 , 5 ) _ J_ 

23,5(1-23,5) '" 4F 2 

Res(M 5 ,z 2 ,4) _ J_ 

22,4(1 - 22,4) " " 4F 2 

„2 



(Sl,4 + «4,5 ' 


-Sl,5j(si,2 +S2,3j 


(Sl,2 + S 2 ,5 - 


Sl,3 

- Sl,5)(«3,4 + S4,5) 


(Sl,4 + S 2 ,5 " 


S3,5 
" «1,5) («2,3 + S3,4) 



«2,4 



+ S 2 ,5 - ^2,4 + s 3,4 - «3,5 



+ Si 4 - Si 3 + S 2 ,3 - S 2 ,4 



+ Si 5 - Si 4 + S 2 ,4 - S 2 ,5 - S2,3 - S3,4 



P- 



1 



[Sl,3 - «1,2 - S 2 ,3 + Sl,5 - Sl,4 + S 2 ,4 ~ S 2 ,5 + S 3 ,5 ~ S3,4 ~ «4,5] • 



(4.57) 



4i?2 4i 72 

Finally we get 
4F 2 M 5 (1) = 

__ (Sj,4 + S 4 , 5 - Si [5 ) (si [2 + S 2 ,3) (Sl,2 + S 2 ,5 ~ Sl,s)(s3,4 + 34,5) ($1,4 + s 2,5 ~ Sl,s) (s 2 ,3 + S 3j4 ) 

Sl,3 S 3j5 

+2S1 5 - Si 2 - Si 4 - S 2 ,3 - S2,5 - S 3 ,4 ~ S 4 ,5- 



S2,4 



(4.58) 



Taking this amplitude on-shell according to (4.53), i.e. setting s± t 5 — > and changing the overall 
sign, we reproduce the parametrization independent physical amplitude (2.34). 

5. More properties of stripped semi-on-shell amplitudes 

The BCFW recursive relations provides us with a Lagrangian-free formulation of the tree-level 
nonlinear SU(N) sigma model in the Cayley parametrization. We can use them similarly as the 
Berends-Giele relations as a tool for the investigation of further interesting features of the stripped 
semi-on-shell amplitudes M2 n +i and J 2 n+i- As we have already mentioned, these features are not 
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universal because of the parameterization dependence of M2 n +\ and J2n+i, however, their implica- 
tions for the fully on shell amplitudes hold universally 11 . In this section we will concentrate on the 
problem of single soft limits (Adler zeroes) and double soft limit of the semi-on-shell amplitudes. 

The presence of Adler zeroes for the on-shell Goldstone boson amplitudes A4 ai '" a2n (pi, . . . ,P2n), 
i.e. validity of the limit 

\lmM aia2 - a2 "( Pl ,p 2 , . • • ,P2n) = 0, (5.1) 

is a well known consequence of the nonlinearly realized chiral symmetry. More generally it is 
an universal (non-perturbative) feature in the theories with spontaneous breakdown of a global 
symmetry. In such theories the amplitudes with one extra Goldstone boson 7r a in the out (or in) 
state vanishes when the Goldstone boson become soft, e.g. 

lim(/ + 7r a (p),out|i,in) = 0, (5.2) 

provided the ir a cannot be emitted from the external lines corresponding to the states |i,in) or 
|/, out). In the SU(N) nonlinear sigma model the Adler zero is present also for the stripped 
on-shell amplitudes M2n(pi,P2, • • • ,P2n) due to the leading N orthogonality relations (2.20) and 
corresponding uniqueness of the decomposition (2.11). However, this property is not guaranteed 
automatically for the semi-on-shell amplitudes M2 n +i and the soft Goldstone boson behavior can 
depend on the parametrization. For instance using the Cayley parametrization, we find for the 
amplitude M% = (p\ ■ ps)/F 2 the Adler zero for soft p\ and p^, however there is no zero for soft P2 
in general when keeping p^ off-shell. For the same amplitude in the exponential parametrization 
(cf. (4.16)) we have no Adler zero at all. As we shall show in this section, for the semi-on-shell 
amplitudes M2 n +i in the Cayley parametrization we can prove, using the BCFW-like relation, the 
Adler zero for half of the momenta (namely for those pj with odd index j). 

The double soft limit of the Goldstone boson on-shell amplitudes .A/pi a 2--- a 2n+2 (j) 1 ^p 2 , • • ■ ,P2n+2) 
is more complicated and has been studied relatively recently in connection with the regularized 
action of the broken generators on the n Goldstone boson states [50]. Motivated by direct inspection 
of the six Goldstone boson amplitude in the nonlinear chiral SU{2) sigma model it was conjectured 
that provided the two soft momenta are sent to zero with the same rate, the following limit holds 

limM abaia2 - a ^(tp,tq,p 1 ,p2,...,P2n) 

= -^E^^ Ca ^'! P I 9 l - Mai ''' a '' 1 " a ' +1,,,a2 "^ 1 ^^---^2n), (5.3) 

2F ~[ Pi- {p + q) 

where / are the structure constants. Analogous statement has been then rigorously proven for 
the tree-level amplitudes in the J\f = 8 supergravity using BCFW relations. In fact, for the on- 
shell amplitudes, the formula (5.3) can be proven non-perturbatively under some assumptions for 
the general enough case of the theory with global symmetry breaking (including the case of chiral 
nonlinear sigma model with general chiral group G) using the symmetry arguments only (cf. the 
PC AC soft-pions theorems [48]) . We postpone the details to the Appendix E . 



Let us remind that the on-shell amplitudes are parametrization independent. 
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In terms of the stripped on-shell amplitudes the relation (5.3) can be rewritten as 

limX 2 n+2(Pl, • • • ,Pi-l,tpi, ■ ■ .,tpj,Pj+l, . ..p2n+2) 



t^o 



1 x { Pi+2 ■ (Pi - Pi+l) Pi-l ■ (Pi - Pi+l)\ . , , . , r . 

;Oj, i+ l -. -, M2n{Pl,---, Pi-l, Pi+2, ■■■P2n+2)- (5.4) 



AF^^\p i+ 2-(pi-p i+1 R-i-(R-ft+i 

In this section we will prove this relation also for the tree-level semi-on-shell amplitudes J2n+i (and 
consequently for M2 n +i) of the SU(N) nonlinear sigma model in the Cayley parametrization using 
suitable form of the generalized BCFW representation. 

5.1 Adler zeroes 

In this subsection we will use the BCFW-like relations (4.52) derived in the previous section and 
prove the presence an Adler zero at M2 n +i when one of the odd momenta, say P21-1, is soft, i.e. we 
will prove that for / = 1, 2, . . . , n + 1 

HmM 2n+ l{pi,P2, ■ ■ ■ ,P2l-2,tp2l-l,P2l+l, ■ ■ ■ ,P2n+l) = 0. (5.5) 

For the fundamental amplitude M 3 (pi,p2,p 3 ) we have explicitly 12 

M 3 (tp 1 ,p 2 ,p 3 ) = M 3 (p 1 ,p 2 , tp 3 ) = —2 t(pi ■ p 3 ) ->• 0. (5.6) 

In the general case the proof of (5.5) is by induction. Let us assume validity of (5.5) for m < n. 
This assumption also means that, taking the cut semi-on-shell amplitude M2 m +i on shell, i.e. for 
p\ 2n+ i — > 0, the Adler zero is in fact present at M2 m +i\ on shell = — ■M.2m+2 for all momenta, i.e. 

]imM2 m+ l(pi,P2,---,tpj,...p2m+l)\on shell = (5.7) 

for all j = 1, ... , 2m + 1 due to the cyclicity of M.2m+2- 

Let us now substitute P21-1 — > tp2i~\ to the right hand side of (4.52). Note that, under such 
substitution, the position of the poles Z2j,2j+2, £2.7-1,2.7+1 and z i ■ become t— dependent. The t— 
dependence of the right hand side of (4.52) is therefore both explicit (due to the explicit dependence 
on P2/-1) and implicit (due to the implicit t— dependence of the poles zp). 

We will now inspect the behavior of the individual terms under the limit t —$■ 0. The first term 
gives finite limit 

(2F 2 ) nP ~ ~* (2F 2 ) nP ^ P21 - 1 ^ ' ^ 5 ' 8 ^ 

As far as the second term is concerned, the individual terms of the sum over j vanish in this 
limit unless j = I — 1. The reason is as follows. For j 7^ I — 1 (the case A in the Figure 9), the 



12 



Note however that for t — > according to (4.18). 



1 . .2 



M 3 (pi,tp 2 ,p 3 ) ->■ ^-fai +P3) 
and therefore the statement analogous to (5.5) for even momenta does not hold. 
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tp2l-l 




P2J-1 



PW+faw) x l?£*3tt2->0 



1 P 2j,2j+2 



P2n+1 \ P2J+3 

Z2j,2j+2P2j+4 



A) j^l-1 




■ P2J-1 



°^^ + tp ^ x ^S,) 2 ^ p^ 2 - ^-^ 



P2n+1 \ P2J+3 

0(t) P2]+i 

B) j =1-1 



Figure 9: Graphical representation of the t —> limit of the second term on the right hand side of (4.52). 
The soft momentum is denoted by dashed line in the case A. In the case B, 0(t) indicates the order of the 
t— dependent 2>y,2j+2- 

kinematical factor p 2 j ■ P2J+2/P2J 27+2 as wen as * ne position of the pole z 2 j >2 j +2 are t— independent 
and because tp 2 i-\ is placed on the odd position in M L h J (^2j,2j+2)) we can safely 13 use the 
induction hypothesis to conclude that 

limM^ 2j+2 \z 2j , 2j+2 )\ P2l _^ = 0. 

For j = 1—1 (the case B in the Figure 9), the kinematical factor p 2 j -p 2 j+ 2 /p2j 21+2 becomes explicitly 
t— dependent and tends to 1/2 for t — >• 0, while M^' J (z 2 j, 2 j+ 2 ) has both explicit (through 
P2j,2j+2 = Z2j,2j+2(P2j + P2J+2) + tP2j+i) and implicit i-dependence. In this case z 2 j t 2j+2 = 0(t), 
as can be seen from (4.43). Therefore, all even momenta in M^ h 3 (z2j,2j+2) are scaled by 0(t) 
factor, in the same way as in (4.18). We can therefore conclude with help of (4.18) that 

limX ia+ ^. W ) p / +2 %y +2 = ^.-i^iU-^o. (5.9) 

The third term on the right hand side of (4.52) can be treated exactly in the same way as the 



^Indeed, in general the momenta pt (z2j,2j+2) and P2j, 2^+2(^2^,2^+2) are t— independent and nonzero. 
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Z2j-l,2j+lP2 

tpu-1 



z 2j-1.2j+lP2j-2 




toj-lM+liz^w+i) X 1 P 2i-l,2i + l-P2i-l,2 i+ l 



P2j-l,2j+l 




A) jVJ.i-1 

0(t) P 2 

0{t)p2j-2 

| tp 2 J-l + 0{t)p2j+P2j+l x 1 P2j-(tP2i-l+P2j+l) . 1 2 I 

l j ' ^ ? (*P2j-l+P2j+P2j+l) 2 (2F2jSP-lp2j-l->0 

P2n+1 " • ' A 0(t)p 2j+ 2 

Pij+3 

b) j = i 

Figure 10: Graphical representation of the t — > limit of the third term on the right hand side of (4.52). 
The soft momentum is denoted by dashed line in the picture A. In the picture B, we show only the j = I 
case, the j = I — 1 case is treated analogously. 0{t) indicates the order of the t— dependent z- 2 j : 2j+2- 

second (see Fig. 10). Also here the individual terms of the sum over j do not contribute with the 
only exception of j = I and j = I — 1 by induction hypothesis applied to M L ' J (z2j-i,2j+i) 
which has for j ^ 1,1 — 1 only explicit t— dependence. In the remaining two cases j = I and j = l — l, 
the explicitly t— dependent kinematical factors V2j-\ 2?+i ' P2j-i 2j+i/p2j-l 2j+i tend again to 1/2 
and within M L ° '' J (z2j-i,2j+i) the even momenta are scaled by Z2j-i,2j+i = 0(t) (see (4.41)) 
and thus (4.18) can be used 14 to conclude that 

Km MW-W+Vf* \ 1 P2j-l,2j + l'P2j-l,2j+l _ x. r s _J 2 | 

\™ Q M L {Z2j - 1 ' 2j+l) ^ 2j _ 12j+1 T 2 " ( j ' 1 + j ' l - l] j2F^ P -^'-^°- 

(5.10) 
The fourth term on the right hand side of (4.52) vanish completely in the limit t — > 0. This is 
easy to see for those terms of the sum over (i,j) for which 15 liiri£_>o z\ ■ 7^ 0. In this case either 
M^ (z i •) or M^ {z i ■) have explicit t— dependence through tp2i—i (which is for M^' {zf-) on 



14 Note that, the odd momenta are t— idependent with the only exception o{p2j-i,2j+l(z2j-i,2j+i)\p 2 4^i-np2jii the 
limit of which is P2j±i- 

15 It is easy to realize that lim t _).o z tj 7^ li m t^o z~ . for generic Pk- 
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odd position) and thus the induction hypothesis in the form (5.5) or (5.7) can be used 16 . By direct 
inspection of (4.40) we find that the only case for which the above argumentation does not apply is 
the case j — i = 4 with i even and i < 21 — 1 < j. Here lim£_>o z~- ^ and so for the "minus" part 
of this (i,j) term we can use the induction hypothesis as above. However, the "plus" part might 
be problematic because 

z+ = - fa'-* • ft"-**;) t + o(t 2 ). (5.ii) 

(P2l-1±2 ■ PTj) 

Using this formula and (4.15) we find after some algebra 

m£ j '>(z+) = M 5 (r«}), • • • toi-i, • • • ,Pi«f )) = 0(t 2 ). (5.12) 

which shows that also the "plus" part has vanishing t — > limit. 

Putting therefore the only nonzero contributions (5.8), (5.9) and (5.10) together we get finally 

limM 2n+ l(pi,P2, • • • ,P2l-2,tp2l-l,P2l+l, ■ ■ ■ ,P2n+l) 
_. / n— 1 n 

which finishes the proof. 

5.2 Double-soft limit 

Let us now study the behavior of the semi-on-shell amplitude J2«+i in the Cayley parametrization 
under the double soft limit, i.e. the case when two external momenta, say pi and pj, are scaled 
according to pij — > tpij and t is sent to zero. In this section we will prove, that for 1 < i < j < 2n+l 

Jim J2n+l(pi, • • • ,P2n+l)|pi-»-tpi,pj-»-tpj 

1 ( (Pi-Pi+2) (Pi-Pi-l) \ T , \fcio\ 

= 6 h i + 1 ^^2 1 i T 1 i T J 2n-l(Pl, • • • ,Pi-l,Pi+2 ■ ■ ■ ,P2n+l)(5.13) 

IF 1 \p i+2 ■ (Pi+l+Pi) Pi-l ■ (Pi+l + Pi) ) 

which has an identical form as (5.4) 17 . The key ingredient of the proof is the generalized form of 
the BCFW representation mentioned in Section 3.2 written for a suitable two-parameter complex 
deformation of the amplitude J2n+i- Such a representation allows us to calculate the double soft 
limit with help of the known behavior of the poles and corresponding residues in this limit. Useful 
information on this behavior can be inferred from the statement (5.5) concerning the Adler zeroes 
proved in the previous subsection. 

The above mentioned deformation of J2n+i can be defined as the following function of two 
complex variables z and t 

Slj{z,t) = J( Pl , . . . ,P2n+l)\p i -^tp i ,p j ->zp j , (5-14) 



16 Let us remind that M ( ^' J \z+ J ) is fully on-shell. 
17 Indeed, 



(Pi ■ Pi+2) (Pi ■ Pi-l) 1 fPi+2 ■ (Pi - Pi+l) Pi-i ■ (Pi - Pi+l 



Pi+2 ■ (Pi+l +Pi) Pi-l ■ (Pi+l +Pi) 2 \Pi+2 ■ {Pi - Pi+l) Pi-l ■ (Pi — Pi+l) 
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therefore 

-^(l, 1) = J2n+l(Pl, • • • ,P2n+l) (5.15) 

Various types of the double soft limit correspond then to various ways of taking the limit (z, t) 



(0,0) in the double complex plane (z,t); the limit (5.13) corresponds to hm t _^o Sf • (t, t) = S*™'- . 



For z — > oo and t > fixed the following asymptotic behavior holds 

S?j(z,t) = O(z ), (5.16) 

as can be easily proved e.g. by induction with help of the Berends-Giele recursive relations (4.7). 
We can therefore write the generalized BCFW relation with one subtraction in the form (3.14) 



Res [SV-j; z k j(t) . 

S?j(z,t) = S?j(a,t) + J2 V m ' j-, . (5.17) 

' J ,J *j-f z-z k ,i{t) z K i{t)-a 

where a ^ z kj i(t) is a priory arbitrary, however, as we shall see in what follows, appropriate choice 
of a can simplify the calculation. 

The poles z k j(t) for k < j < I correspond to the conditions p\ i\pi-ttpi# —*zp ■ = 0) or explicitly 

z k ,m = - f u ^f "° . (5.i8) 

2(Pj • Pk,l)\pi->tpi 

The residues at the poles z k i(t) factorize 

Res (Sfj; z k ,i(i)) = —. r, [J2n+l~(l-k)(Pl, ■ ■ ■ ,Pk-i,Pk,l,Pl+i, ■ ■ ■ ,P2N+l) 

z \Pj " Pk,l)\pi->-tpi 

xM t _ k+1 (p k , . . . ,p0lpi-»-tpi^-»-2pJU-»-zfe,j(t)> (5-19) 

where Mi_ k+ i is the cut amplitude (4.30). Namely the latter two formulae along with (5.5) contain 
sufficient amount of information for the calculation of the double soft limit. 

Let us first assume i < j where i is odd and j arbitrary. This choice is a technical one, and as we 
shall see, the general case can be easily obtained using the symmetry properties of the amplitude. 
In what follows we set a = 1 in (5.17), the double soft limit then simplifies to 



Res SP,-;z fcl j(t t _ 1 

S?f = lim Sf , (t,t) = lim V V , . ' — , (5.20) 

l ' J t^^ K '> t^o^ t-z k ,i(t) z k>l (t)-V 

where we have used the existence of the Adler zero for S™Al,t) = J2n+i(pi, ■ ■ ■ , tpi, . . . ,P2n+i) and 
* odd (cf. (5.5)). 

For generic p r there exist a finite limit 

z k>l (0) = \imz k , l (t)^l (5.21) 

t-¥0 
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In fact the only nonzero contributions to the right hand side of (5.20) stem from the cases for which 
Zk,l(0) = 0. Indeed, for Zk,i(0) /Owe get for the corresponding contribution 

-lim Res (S^.; %/(*)), (5.22) 



**,j(o)CM°) 

and, according to (5.5), on the right hand side of (5.19) we get either 



lim[Mi- k+1 (p k , • • • ,Pi)lpi-+^Nj»j-^«p J -]U-+« fc ,j(t) = ° ( 5 - 23 ) 



for k < i < j < I or 



lim J 2n+ i-(i- k )iPh ■ ■ -,tPi, ■ ■ ■ ,P(k, l)(t),pk+i, ■ ■ ■ ,P2n+i) = (5.24) 

for i < k < j <l. In both cases the complementary factor has finite limit and therefore 

limRes(5^;z M (i))=0. (5.25) 

Let us therefore discuss the contributions form the poles for which Zk 7 i(0) = 0. Note that, for 
generic p r such a pole does not exist provided j > i + 2. We can therefore immediately conclude 

Slf = for j > i + 2. (5.26) 

What remains are the following two alternatives for which the three-particle poles Zk t i(t) with 
I = k + 2 can vanish in the limit t — > (see Fig. 11) 

1. j = i + 1 and either k = i or k = i — 1. In this case either 

Pi-l,i+ilpi->*Pi,Pj— >o — >-Pi_i = (5-27) 

or 

Pi,i+2\pi^t Pi ,p :j ^o ->■ Pi+2 = (5.28) 

2. j = i + 2 and fe = i, in this case 

Pti+2\pi^t Pi , Pj ^o = pf+i = 0. (5.29) 

In what follows we will discuss separately the cases j = i + 1 and j = i + 2. Let us first study the 
double soft limit of two adjacent momenta, i.e. j = i + 1 where i is odd. We will investigate the 
contributions of individual poles Zk,i(t) on the right hand side of (5.20) separately. In this case we 
get for i > 1 only two potentially nonzero contributions (i.e. (5.28) and (5.27)) to the right hand 
side of (5.20), namely 

Res fem;^-M+i(0) t-i . Res fe+i;^+2(*)) t-i 



S™},-, = lim — : — : — f- + lim ■ 

' + *^0 t - Zi-i :i+1 {t) Zi- lli+1 (t) - 1 *^o t-z iji+2 (t) Zi, i+2 {t)-l 



(5.30) 
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#<-l,i+l(t) — ► 




2n+l 



i-2 



Pi-i + tpi + Zi_ u+ i{t)p i+1 



i + 2 



i-2 




i- 1 



\ i + 1 



(PrK-l) 
Pi-r(p»+K+l) 




2-1 



2?7 + 1 ' • * * + 2 



-K*,*+-2(*) 




i- 1 



Pi+2+«Pi +\i+2{t)Pi+l 



i-1 



(M;, )- * « + 1 _^. (Pi-Pi+2) 




K+2'fe+Pi+l) 



j + 2 



2n+l * • * * + 3 




J + 2 



2n + l • • 



i + 3 



^, i+2 (i) 




i- 1 



Pi+l+*Pi+Z,,i+2(*)Pi+2 



(M 3 u+ 
\ 



i + l 

+ 2 



2n+l 



i + 3 



Figure 11: Graphical representation of the t — >• limit of the three cases (5.27), (5.28) and (5.29) for which 
Zk,i(t) — > 0. The soft momenta are denoted by dotted lines. The multiplicative factors K k ,i(t) stays for 
(t'-l)/(t-z kJ (t))(z k ,i(t)-l). 



We get for the poles ^»-i,i+i(t) and 2^1+2 (£) 

^fc,fc+2(*) 



Pfc,fc+2lpi^*Pi,Pj^0 (p, • p r ) /-y^2\ 



2(pj • Pfc,fc+2) |pi->tpi (Pj • Pr) 

where either r = i + 2 (for k = i) or r = i — 1 (for A; = i — 1), and as a consequence, 

1 i-1 _ 1 (pj -p r ) 



(5.31) 



We have further 



* - z kl k+2(t) z k)k+2 {t) - 1 £ p r • (pj + p. 

Pk,k+2{t) = tpi + Z Kk+2 {t)p,j +p r -+p r ^ 



(i + o(t)). 



(5.32) 



(5.33) 



and therefore in both cases 



lim J 2n -l(Pl, • • • ,Pfc-l,Pfe,ifc+2(*),Pfe+3, ■ ■ ■ ,P2ra+l) = ^2n-l(Pl, • • • ,Pi-2,Pi-l,Pi+2, • • • ,P2n+l)- 



(5.34) 
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For the remaining ingredients of the formula (5.19) we get 

M 3 (tpi,Zi : i +2 (t)pi + i,pi + 2) = -0t(pi ■ p i+2 ) (5.35) 

Mz{pi_i,tpi, Zi-i t i + i{t)p i+ i) = —^Zi-. 1>i+1 (t){pi-i -Pi+i) = -t—^ipi -pi-i)(l + 0(t)). (5.36) 
Inserting this into the formulae (5.19) and (5.30) get finally for i > 1 



S i,i+l = 7T^> ( ~ 71 T-7T ~ Z — T-7T ) J 2n-l{Pl, ■ ■ ■ ,Pi-2,Pi-l,Pi+2, • • • ,P2n+l)- 



,„.{) _J_ ( (Pi-Pi+2) (Pi -Pi-l) 

2F 2 \pi+2 ■ (Pi+i +Pi) Pi-i ■ (Pi+i +Pi) 

(5.37) 

In the same way, for i = 1 only the first term on the right hand side of (5.37) contributes. 

Let us proceed to the case 2. when j = i + 2 and Zi^+iit) — > for t — > is the only pole which 

can give nonzero contribution to (5.20). In this case we have 

l ' t+2 t-*o t-z hl+2 (t) Zi, i+2 (t)-l v ; 

The formulae (5.31, 5.32, 5.33, 5.34) are still valid with r = i + 1, but now we have 

M3(tpi,p i+ l, Zi,i+2(t)pi+2) = -T^tZ iyi+ 2(t)(pi ■ Pi+2) = 0{t 2 ). (5.39) 

infl 



which implies S i \ +2 = 0. 

To summarize, we have for k > 



lim J 2n +l(Pl, ■ ■ ■ ,P2k,tp 2 k+l, ...,tpj,... ,p 2n +l) 



1 T , s ( {P2k+1 ■ P2k+3) (P2k+1-P2k) 



^j,2k+2TT^2^n-l{pi, ■ ■ ■ ,P2k,P2k+3, ■ ■ ■ ,P2n+l) 



Zr \P2k+3 • [P2k+2 + P2k+l) P2k • {P2k+2 + P2k+1 

(5.40) 
and for k = 

lim J 2n+ i(tpi, ...,tpj,... ,p 2n+ i) = 5j,27TEvl 7T7 ^2n-l(P3, • • • ,P2n+l)- 

t^o 2F 2 (p 2 ■ p 3 ) + (pi • p 3 ) 

As it is clear from the above discussion, the "asymmetry" of the latter result stems from the fact 
that P2n+2 is off-shell and therefore the three-particle pole corresponding to (p 3 +P4 + . ■ ■+P2n+\) 2 = 
{P2n+2 — Pi — P2) 2 — > P 2n +2 7^ does not contribute. 
Because 

J(l, 2, . . . , 2n + 1) = J(2n + 1, 2n, . . . , 2, 1), (5.41) 

we get for j < 2k + 1 

J2n+l(Pl,--- ,tpj ... ,P2k,tp2k+l, ■ ■ ■ , ,P2n+l) = ^2n+l(P2n+l, ■ ■ ■ ,tP2k+l,P2k, ■ ■ ■ ,tPj,--- ,Pl)- 

(5.42) 
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On the right hand side of this identity the momentum P2k+i stays on the odd position and thus 



lim J2n+l(pi, ■■-,tpj---,P2k, tp2k+l, ■ ■ ■ ,P2n+l) 



x 1 T , s( (P2k-P2k-l) (P2k-P2k+2 

0j,2k7TB2 Jl 2n-l{Pl ■ ■ ■ ,P2k-l,P2k+2, ■ ■ ■ ,P2n+l) I 1 : T + 



2-F 2 \ P2k-1 ■ [P2k + P2k+l) P2k+2 ■ \P2k + P2fc+1 

(5.43) 
Putting (5.40) and (5.43) together the final result (5.13) follows. 

6. Summary and conclusion 

We have studied various aspects of the SU(N) chiral nonlinear sigma model which describes the 
low-energy dynamics of the Goldstone bosons corresponding to the spontaneous chiral symmetry 
breaking SU(N) x SU(N) — > SU(N). As we have shown, the tree-level scattering amplitudes of 
the Goldstone bosons can be constructed from the stripped amplitudes, which are identical as those 
of the U(N) chiral nonlinear sigma model. It is therefore possible to use this correspondence and 
to investigate both the SU(N) and U(N) cases on the same footing. Especially we are allowed 
to choose any parametrization (field redefinition) of the chiral unitary matrix U(x) entering the 
Lagrangian from the wide class of parametrizations admissible for the extended U(N) case, because 
the fully on-shell stripped amplitudes do not depend on the parametrization. For the direct calcu- 
lation of the flavor ordered Feynman graphs, the most convenient choice proved to be the minimal 
parametrization (2.31), which we have chosen in order to calculate the on-shell amplitudes up to 
10 Goldstone bosons. 

The proliferation of the Feynman graphs with increasing number of the Goldstone bosons 
call for alternative methods of calculation. The more efficient method is based on the Berends- 
Giele recursive relations for the semi-on-shell amplitudes, but due to the infinite number of the 
interaction vertices in the Lagrangian of the nonlinear sigma model, the number of terms necessary 
to evaluate the n— point amplitude grows much faster (exponentially) with n than for the case of 
the power-counting renormalizable theories (where the growth is polynomial). 

The BCFW recursive relations could make the calculation of the on-shell stripped amplitude as 
effective as for the renormalizable theories at least as far as the number of terms (which is in both 
cases related to the number of factorization channels) is concerned. However, the standard way of 
the BCFW reconstruction is not directly applicable for the nonlinear sigma model because of the bad 
behavior of the BCFW deformed amplitudes at infinity. We have therefore proposed an alternative 
deformation of the semi-on-shell amplitudes based on the scaling of all odd or all even momenta, for 
which we were able to prove exact results concerning the behavior of the semi-on-shell amplitudes 
when the scaling parameter tended to zero. Using the Berends-Giele recursive relations we were 
able to prove this scaling properties for general n— point amplitude. An essential ingredient of the 
proof was the fact that the semi-on-shell amplitudes (unlike the on-shell ones) are parametrization 
dependent and we could therefore make an appropriate choice of the parametrization (the Cayley 
one). We have then used these exact scaling properties for a generalized BCFW reconstruction 
formula (with one subtraction) which determines fully all the semi-on-shell amplitudes in the Cayley 
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parametrization including the basic four-point one. Putting then the semi-on-shell amplitudes on- 
shell we reconstruct simply the parametrization independent on-shell amplitudes. In contrast to 
the standard BCFW relations our procedure is not restricted to d > 4 space-time dimensions. 

The BCFW recursive relation are also a suitable tool for investigation of the properties of the 
amplitudes. We have illustrated this in two cases, namely we have proved the presence of the Adler 
zero and established the general form of the double soft limit for the semi-on-shell amplitudes in 
the Cayley parametrization. 

The existence of BCFW recursion relations for power-counting non-renormalizable effective 
theory as the SU(N) chiral nonlinear sigma model gives an evidence that the on-shell methods can 
be used for much larger classes of theories than has been considered so far. It also indicates that the 
SU(N) chiral nonlinear sigma model is rather special and deeper understanding of all its properties 
is desirable. For future directions, it would be interesting to see whether the construction can be 
re-formulated purely in terms of on-shell scattering amplitudes not using the semi-on-shell ones. 
Next possibility is to focus on loop amplitudes. As was shown in [20] the loop integrand can be 
also in certain cases constructed using BCFW recursion relations, it would be spectacular if the 
similar construction can be applied for effective field theories. 
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A. General parametrization 

In this Appendix we will discuss a very general class of parameterizations of the U(N) sigma model 
originally studied in [1], which is suited for a derivation of the stripped Feynman rules. Within this 
class the field U{x) £ U(N) is expressed in the form 

oo , . \k 

U = Y,^k{V2^A (A.l) 

where (j) = t a (j) a , <j) a are the Goldstone boson fields, t a are the U(N) generators normalized according 
to (t a t ) = 5 and a k are real coefficients. These coefficients are not completely arbitrary, because 
the unitarity condition U + U = 1 implies the following constraint 

n 

^2a k a n - k (-l) k = S nfi . (A. 2) 

fc=o 

For n = we get Oq = 1 and without lose of generality we can set ao = 1 . In order to preserve the 
correct normalization of the kinetic term and to keep the interpretation of F as the decay constant 
for the fields (ft a we have to fix also a± = 1. 
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For n odd the relations (A. 2) are satisfied automatically while for n = 2k we can solve them 
for ci2k and get a recurrent formula for the even coefficients expressed in terms of the odd ones 

I _-\\k k~l 

«2fc = 5 °fc ~~ ^(-l) i Oj02fc-i- ( A - 3 ) 



This gives up to k = 3 



1 2 1 

ao = — Oi = — 

2 x 2 

1 2 1 

a 4 = --a 2 + aia 3 = -- + a 3 

z o 



1 2 111 

-a 3 + aia 5 - a 2 a 4 = — - -a 3 + -< 
z lb z z 



a 6 = 7T a 3 + a l°5 — a 2«4 = — — -^3 + 7T«3 + «5 (A. 4) 



The explicit solution of the recurrent relations (A. 3) to all orders can be easily found by means 
of the following trick. Let us introduce the generating function f{x) of the above coefficients a/j 

oc 

f(x) = Y,akX k . (A.5) 

fc=0 

The relations of unitarity with the initial conditions oo = 01 = 1 are then equivalent to 

f(-x)f(x) = 1, /(0) = 1, /'(0) = 1 (A.6) 

which represents a functional equations for the generating functions f(x). Let us define f±(x) to 
be the even and odd part of f(x), i.e. f±(x) = (f(x) ± /(— x)) /2. From (A.6) we get then 

f + {x) 2 -f-{x) 2 = l (A.7) 

or finally 

f + (x) = V1 + /-W 2 - (A.8) 

The formal series expansion of both sides of the last equation at x = gives the solution of the 
recurrent relations (A. 3), i.e. the explicit expressions for a 2 fc in terms of an infinite number of free 
parameters fl2fc+i- The general solution of the functional equation (A.6) is then 

f(x) = /_ (x) + V1 + /-M 2 (A.9) 

where f-{x) is arbitrary odd real function analytic for x = satisfying /'(0) = 1. The minimal 
parameter-free solution corresponds to the choice a2k+i = for k > 0, i.e. /™ m (x) = x and 



/min(z) = X + \A + X 2 (A. 10) 

i.e. for fc > 1 



i v x 



# = ( „!,?,! C fc -i, (A.ll) 
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where 




0.-^1" I 



Jexp (X J 


= 


c j; 




/Caylcy(x) 


= 


1 + (x/2) 
1 - (x/2) 


exp 
«fc = 


1 






Cayley _ 
a k 




1 


1 

r>fc-1 



are the Catalan numbers. 

Another frequently used choices are the exponential and Cayley parameterizations correspond- 
ing to /cx P (a;) and /cayicy(x) respectively, where 

(A.13) 
(A.14) 

1 - \X/ Z) 

or in terms of the coefficients a k 

(A.15) 

(A.16) 
J- -r uk,o *~ ' 

These two parameterizations can be understood as minimal parameter-free variants with respect 
to other two possible forms of the general solutions of the functional equation (A. 6), namely 

f(x)=expg(x) (A.17) 

and 

«*> = Sf§) (A - 18) 

where g(x) and h(x) are arbitrary real functions analytic for x = for which 

g(x) = -g(-x), (A.19) 

g(0) = 0, g'(0) = 1 (A.20) 

and 

h'(0) = -h(0) ± 0. (A.21) 

As was proved in [1], for N > 2 the only parametrization from the class (A.l) admissible also 
for SU(N) sigma model is the exponential one. The reason is that, under the general axial SU(N) 
transformation 

OO f . \ k CO f . \ k 

U(xy = ^2a k (V2^') =U A J2ak[V2j<t>) U A (A.22) 

fc=o ^ ' fc=0 ^ ' 

which defines corresponding nonlinear transformation of the matrix of the Goldstone boson fields 
<j) = ^2 a= i <j) a t a the SU(N) condition for the trace ((f)') = is not preserved unless a^ = l/k\. 
Of course, in the case iV > 2 we can use different admissible parameterizations of SU(N) which, 
however, do not belong to the class (A.l) (see e.g. [47]). 

Let us now find the stripped Feynman rules. Using the general parametrization (A.l) we can 
write the Lagrangian of the nonlinear U(N) sigma model in the expanded form 
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p2 °° 

£( 2 ) = —(dU ■ dU + ) = Y, v n ,m{ d ^ n ■ W">- (A.23) 

n,m=0 

where we get for w„ irn after some algebra (and using the unitarity condition (A. 2)) 

v n , m = (1 + (-l) n+m )^4+^T E «ka m+n+ 2-k{-l) k+ \k - 1 - m) (A.24) 

A;=0 

Therefore only the terms with even number of fields survive, explicitly 



£ (2) = E45 +2 (A-25) 



t ,(2) 

ra=0 

where 

2n 

45+2 = E «*,2n-fc<0#* • 300 2n - fc ) (A.26) 

fe=0 

The usual Feynman rules for the vertices can be easily obtained as a sum over permutations 

C;2 ,a2 " +2 (P1>P2, • • • ,P2n + i;P2n+2) = "2" +1 £ <i<MD . . . tM***)) 

2n 

X ^ Vfc,2n-fc(Pa(l) • P<r(l)+fc+l) (A.27) 

fc=0 

The stripped Feynman rule then follows in the form 

2n 2n+2 
V 2n +2{Pl,P2, ■ ■ ■ ,P2n+i;P2n+2) = -2™ +1 ^ ^ V kj2n -k{Pi ' Pi+k+l) (A. 28) 

fc=0 i=l 

Inserting (A. 15) into (A.24) we get after some algebra for the exponential parametrization 

ex P _ (-!)" (~l) fc /^2iA 
v fe,2n-fc " 2F2« (2n + 2)! ^ fc J ' ( 9) 

while for the Cayley parametrization we have v 2k+i Y 2n-2k-i = ^ an( ^ 

Cayley _ l" 1 )" * /a on\ 

u 2k,2n-2k - 2i?2n 2 2n+l ' t^.OUj 

Similar calculations can be made also for the minimal parametrization, but the result is much more 
lengthy and we will not need it explicitly Instead we will rewrite the Feynman rules for the vertex 
V2n+2 with 2n + 2 external legs in terms of the variables 

2 



Sij =p {J (A. 31) 
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where 1 < i < j < 2n + 1 and 

3 

Pi,j = ^Pk (A.32) 

k=i 

Here we identify 

S2n+2,2n+2+k = s k+l,2n+l (A. 33) 

Si,2n+2+k = Sfc+l,i-l- (A. 34) 

The scalar products (pi ■ pj) can be then expressed as 

(Pi ■ Pi) = Si.i (A.35) 

(pi ■ Pi+i) = p(sj,i+i - Sj,i - Sj+i^+i) (A.36) 

and for k > 2 

(Pi " Pi+fc) = „(Si,i+fc _ Si,i+fe-l + Sj+l,i+fc-l - «i+l,i+fc)- (A. 37) 

On-shell we get s$ $ = and si 2n+i = 0- The stripped Feynman rule in these variables can be 
written in the form valid for n > 1 

/ r, \ n n 2n+2 

V 2 n+2{Si,j) = (-1)" ( -^2 J 5Z^fc,n J^ ^i+fc (A.38) 

^ ' fc=0 i=l 

where 

™0,n = ("l) n 2F 2n (2t> ,2n " Ul,2n-l) (A.39) 

Wfc,n = (~l) n 2F 2n (2v k ,2n~k - Vk-l,2n+l-k - Vk+l,2n-l-k) for k < Tl (A.40) 

^n,n = (-l) n 2F 2 "K, n - t; n _i in+ i). (A.41) 

Within the general parametrization we get from (A. 24) and (A. 2) after some algebra 

(-l) k 
w k ,n = , : x a k+1 a 2n +i-k- (A. 42) 

1 + Ofcn 

For the above special cases this reads for N > 1 



w ex P (-l) fc 1 /2n + 2\ 

U,fc ' n _ 1 + **« (2n + 2)! ^fe + lj ( j 

Caylcy _ (- 1 ) l _ /a ^ 



<n = <k m n = (A.45) 

«>&,„ = T— ^ ^C fc C n _ fc _L (A.46) 

l + <52fc+l.n 2 /n 



>2k+l,n 

-3 



Note that, for the minimal parametrization the coefficients w™^ 1 at s^i = pj vanish, therefore the 
stripped Feynman rules for vertices do not depend on the off-shellness of the momenta in this case. 
This fact has been observed already in [11] without calculating the explicit Feynman rules. 
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B. More examples of amplitudes 

The eight-point amplitude is 

8F 6 .M(1,2,3,4,5,6,7,8) = 

_ 1 (Sl,2 + S2,3)(S1,4 + 54,7)(S5,6 + 86,7) (^1,2 + gg^gKfM + S^ 5 )(s 6 j + 5 7 ,; 

2 Si 385 7 Si 3 S 6 ,8 



51,3*5,7 

(Sl,2 + ■S2,3)(S4,5 + 84,7 + 85,6 + 85,8 + 86,7 + 87,8) 
81,3 



+ 2ai >2 + 2«i,4 + cycl (B.l) 



and graphically in Fig. 12. Finally the ten-point amplitude is given by 




Figure 12: Graphical representation of the 8-point amplitude (B.l) with cycling tacitly assumed. 



16F 8 M(1, 2, 3, 4, 5, 6, 7, 8, 9, 10) 



Sl,2 + 82,3 
s l,3 



1 (Si, 4 + S 4i 9)(s 5i 8 + S6,9)(S6,7 + S 7) g) . 1 («1,4 + 84,5)(si,8 + S^ 9 )(s e j + S 7 , 8 ) 



S5,9S6,8 



+ 



Sl,5S6,l 



_ 1 (Sl,8 + 84,9X84,5 + S5,8)(*6,7 + g^g) (si,4 + 84,5X81,6 + S 6 , 7 )(si,8 + 88,9) 
2 8 4 , 8 S 6 ,8 81,551,7 

(51,4 + S4,5)(S1,6 + S 6 , 9 )(87,8 + 8 8 ,g) (81,8 + S4,g)(s4,7 + S 5 ^)(s 5 ,6 + S 6 , 7 ) 



+ 



Sl,5S7,9 84,885,7 

(5i, 6 + 54,9)(84,5 + 8 5 ,6)(87,8 + 8 8 ,g) _ 1 (81,4 + 5i, 8 + 84,5 + 5 4 , 9 + 5 5 , 8 + 5 6 , 9 )(s 6 ,7 + 5 7 , 8 ) 

2 



84,687,9 

(5i,8 + 54, 9 )(5 4 , 5 + 54,7 + 5 5 , 6 + 5 5 , 8 + 5 6 ,7 + 5 7 , 8 ) 

84,8 
(81,4 + 8l, 6 + 54,5 + S 4 , 7 + S 5 , 6 + S6,7)(8l,8 + 88,9) 

81,7 
(51,4 + 5i, 6 + 84,5 + 84,9 + 85,6 + 86,9) (87,8 + 88,9) 

87,9 
(51,4 + 5 4 , 5 )(5i, 6 + 5i,8 + 86,7 + 86,9 + 87,8 + 88,9) 

81,5 
(51,4 + 5 4 , 9 )(s 5 , 6 + 5 5 , 8 + 5 6 , 7 + 5 6 , 9 + 57,8 + 5 8 , 9 ) 

85,9 



86,8 



+ 2si, 4 + si, 6 + 2si, 8 + 2s 4 , 5 + S 4 , 7 + 2s 4 , 9 + 2s 5 , 6 + 85,8 + 2s 6 , 7 + 8 6 , 9 + 2s 7 , 8 + 2s 8 , 9 | 
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Figure 13: Graphical representation of the 10-point amplitude (B.2) with cycling tacitly assumed. 



1 (si,2 + 51,4 + S2,3 + ^2,5 + ^3,4 + S4, 5 )(si [6 + $1,8 + ^6,7 + ^6,9 + ^7,8 + 38,9) 

2 Sl,5 

+ 5si j2 + 2si ( 4 + cycl 
with one-to-one correspondence with Fig. 13 



(B.2) 



C. Relative efficiency of Feynman diagrams and Berends-Giele relations 

In this appendix we review the solution of several types of recursive relations which count the 
number of ordered Feynman graphs needed for the semi-on-shell amplitude .7(1,2, ... , n) in the 
nonlinear sigma model and related toy models. 

C.l Number of the Feynman graphs 

Let us start with the case of nonlinear sigma model, i.e. with the case with infinite number of 
vertices in the interaction Lagrangian. The above recursive relations, which determine the number 
/(2n + 1) of the (flavor ordered) Feynman graphs which contribute to J(l, 2, . . . , 2ra + 1), are tightly 
related to the Berends-Giele relations (4.7). Indeed, after making the following substitution to (4.7) 



J(l,2,...,2n + l)-)-/(2n + l), 



£*2n+2 



-> 1, i^2fc+i -> 0, iV 5 



2fc+2 



(C.l) 
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the individual terms on the right hand side just count the number of Feynman graphs generated 
from these terms by the iterations of the recursive procedure. As a result we get for f(2n + 1) the 
following recursive relation 

n 2k+l 

/(2n + l) = EE X\f(?n l + l), (C.2) 

k=\ {m} i=i 

with the initial condition /(l) = 1. In the above formula the sum over {n{\ is constrained by the 
requirement 

2fc+l 2fc+l 

E (2m + 1) = 2n + 1 ^ J2m = n-k (C.3) 

i=l i=l 

i.e. it corresponds to the sum over all possible decompositions of ordered set of 2n + 1 momenta to 
non-empty clusters with odd number of momenta in each cluster (cf. (4.8) and Fig. 3), i.e. more 
explicitly 

n 2fc+l 

/(2n + l) = E E X\f{1n l + l), /(1) = 1. (C.4) 

k=lJ2i n i= n —k * = 1 
Standard method for solution of this type of recursive relation is based on the generating function 
defined as 

oo 

A(x) = yf(2n + l)x n . (C.5) 

n=0 

The recursive formula (C.4) implies the following equation for A(x) 

00 A 3 

A = 1 + y x k A 2k+l = 1 + -^— i (C.6) 

^— ' 1 — xA z 

k=l 

° r BE 

x = 5 = "T^T (C7) 

(B + lf{2B + l) g(B) 

where B = A — 1 and g(z) = (z + l) 2 (2z + 1). In this form, the problem is prepared for the 
application of the Lagrange-Biirmann inversion formula 

00 T n An— I °° T n An— 1 

*(*) = E ^i^i9(zru =0 = E ^ a^ ( * + l ^ 2z + l T\*=«- ( c - 8 ) 

n=0 rz=l 



After straightforward algebra with help of Leibnitz rule we get for n > 1 

/(2 „ + 1) = H!lg( » )W 2 -> = ir-Vifi-».-^i), (C.9) 



n ^U+l U V 2 

fc=0 



where 2Fi(ct-/3,y,z) is the hypergeometric function. In the same way one can solve the recurrence 
relations for the number of ordered Feynman graphs for the semi-on-shell amplitudes J(l, 2, . . . , n) 
in the cases when only quadrilinear vertices ("0 4 theory"), only trilinear vertices ("i^ 3 theory") or 
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n 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


h(n) 


1 


2 


5 


14 


42 


132 


429 


1 430 


4 862 


16 796 


h+i{n) 


1 


3 


10 


38 


154 


654 


2 871 


12 925 


59 345 


276 835 


A(2n + 1) 


3 


12 


55 


273 


1 428 


7 752 


43 263 


246 675 


1 430 715 


8 414 640 


/(2n+l) 


4 


21 


126 


818 


5 594 


39 693 


289 510 


2 157 150 


16 348 960 


125 642 146 



Table 2: Number of flavor ordered Feynman graphs for J(l, . . . , n) and J(l, . . . , 2n + 1) in the models of 
the type (f> 3 , (f> 3 + </> 4 , 4> A and nonlinear sigma model. 



both trilinear and quadrilinear vertices ("0 3 + 4 theory") are present in the Lagrangian. In the first 
case, similarly to the nonlinear sigma model, only J(l, 2, . . . , n) with n odd can be nonzero, while 
in the remaining two cases J(l, 2, . . . , n) both parities of n are generally allowed. Let us denote the 
number of the Feynman graphs for J(l, 2, . . . , n) as /4(n), f2(n) and /3+4(n) respectively. We get 
the following recurrence relations 



/ 4 (2n + l) 



Y, /4(2m + l)/ 4 (2n 2 + l)/ 4 (2n 3 + 1) 

ni+n2+ii3=n— 1, rii>0 

hi n ) = E h{ni)h{n2) 

ni+ri2=n, «;>! 

E h+4{ni)f3+4(n 2 ) 



/3+4(n 



ni+n2=n, n;>l 

+ £ 



./ 3+4 ("1 ) /3+4 (n 2 ) /3+4 (n 3 ) 



ni+n2+n s =n, n;>l 

with initial conditions /j(l) = 1, j = 3, 4, 3 + 4. The corresponding generating functions 

oo oo 

M{x) = E ^ 2n + ^^ ^3,3+4(x) = J] /3,3+4(n)x" 

n=l 



n=0 



then satisfy 



A 4 = l + x^, ^ 3 = x + ^l, A 3+4 = x + A^ d + A* 



1 3+4 "i" ^3+4- 



In the second case we get 

A 3 (x) 
and therefore 



1 - VI - 4x 1 



i-E 



1/2 



(-4x)* 



n=0 



n \ n — 1 



^n-l 



where C n are the Catalan numbers. In the first case, writing 

_ A A - 1 _ B 4 

x ~ TT " (Ba + 1) 3 



(CIO) 
(C.ll) 

(C.12) 
(C.13) 



(C.14) 



(C.15) 



(C.16) 



(C.17) 



(C.18) 
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and using the Lagrange-Biirmann inversion formula we get for n > 

1 H n_1 1 / S?7 \ 

« 2 " + 1 » = S^ (2 + »» fc l- = STT ( „ J ■ < C19) 

In the third case, we get from 

x = ^3+4 (1 - A 3+4 - AQ (C.20) 

and using the Lagrange-Biirmann inversion formula 

f , n) = i d"- 1 / i y = (-i)- d*- 1 / i y / i y 

n! dz n_1 \1 — z — z 2 J n! dz ra_1 \£i — z/ \^2 — -2/ 

(where z\ = —<fi, Z2 = <p ~ 1 an d ^ = (1 + v5)/2 is the Golden ratio) the result 

w^(- 1 r^g("-: +fc )( 2( ";_V*)(i^' 

4 \«-i / i 



I (n--J 2 Fifl-n J n,2-2n; T ^l. (C.22) 

The first twelve members of the above sequences are illustrated in the Table 2. 

C.2 Efficiency of the Berends-Giele relations 

We can compare this with the number of terms generated by Berends-Giele recursion. For the 
nonlinear sigma model, the number of terms on the right hand side of (4.7) is just 

n n / . \ 

*(2n+i)=x;E i =Eu-fe )= f ™- 1 ^ 

fc=l {nj fc=l V / 

where 

F n = ^=( ( p n -(4>-m (C24) 

are the Fibonacci numbers and <f> = (1 + \/5)/2 is the Golden ratio. Therefore, using the known 
results for J(l, 2, . . . , 2m + 1) with m < n at each step, we need to evaluate altogether 

b(2n + 1) = ± t(2m + l) = j= U^ ~ (<t> ~ ^TT^y) ~ U ^ 

m=\ V \ V W I / 

terms in order to calculate J(l,2, ...,2n + 1) using the Berends-Giele recursion. We show the 
sequences t(2n + 1) and b(2n + 1) in the first and second row of Tab.l respectively 

In the same way we can calculate analogous numbers tj(n) and bj(n) for j = 3, 4, 3 + 4, i.e. for 
"</> 3 theory" , "0 3 theory" or "</> 3 + </> 4 theory". For instance, for £4(2?! + 1) we have (see Tab. 1 for 
numerical values) 



£ 4 (2n + 1) = I U + l j , 6 4 (2n + 1) = £ £ 4 (2m + 1) = ^ 

\ / m=l 



n(n + l)(n + 2) (C.26) 

6 

Note the exponential growth of t(2n + 1) and 6(2n + 1) with increasing n in contrast to the only 
polynomial growth of £ 4 (2n + 1) and 6 4 (2ra + 1) . 
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D. Other example of scaling properties of the semi-on-shell amplitudes 

In this appendix we prove the following scaling limit 

lim J2n+l(tPl,P2, tp3,P4, ■ ■ ■ , tp2r-l,tp2r,tp2r+l, ■ ■ ■ ,P2N,tp 2n +l) = (D.l) 

which is valid for for n > 1. Let us note, however, that 

J3(tpi,tp 2 ,tp 3 ) = J3(pi,P2,P3) ¥= 0. (D.2) 

On the other hand, for N = 2 we get by direct calculation 

lim J 5 (tpi,tp 2 ,tp3,p4,tp 5 ) = lim J 5 (tpi,p2,tp 3 ,tp4,tp 5 ) = (D.3) 

t— !>0 t-tO 

and we can therefore proceed by induction based on Berends-Giele relations almost exactly as in 
the case of the proof of (4.17). The only modification here is that, along with the "dangerous" 
contributions without blocks J(jk + 1, . . . ,jk+i) where jk is even and jk+i — jk > 1 attached to 
the odd line of the vertex V m +\ (provided at least one such a block is present, the contribution 
vanish either by the induction hypothesis or by (4.17) ) we have to discuss separately new type 
of "dangerous" terms with building block J(p2r-i,P2r,P2r+i) (this block does not vanish due to 
(D.2)). The "old" dangerous terms do not in fact contribute as was already discussed within the 
proof of (4.17). The "new" dangerous terms have the following general form form 

— lV 2 k+2(Pl,P2,2j l ,P2j 1 +l, ■ ■ ■ ,P2Ji+2,2r-2,P2r-l,2r+l,P2r+2,2j l+1 ■ ■ ■ 

P2N+2 

■ ■ ■ ,P2j k _ 1 ,2n,P2n+l, ~Pl,2n+l) 

X J(pi)J(2, . . . , 2j 1 )J( P2n+1 ) . . . J{2j l + 2, . . . , 2r - 2) 

Xj(P2r-l,P2r,P2r+l)J(2r + 2, . . . , 2j i+1 ) ■ ■ ■ J(2j fe _i, . . . , 2n)J(p 2n+ i). (D.4) 

Note that, P2r-i,2r+i is attached to the odd line of the vertex V2k+2 and scales as 

P2r- l,2r+l — > tp2r-l,2r+l (D-5) 

i.e. in the same way as the remaining momenta attached to the odd lines of the vertex. The vertex 
being proportional the squared sum of the odd line momenta scales therefore as 0(t 2 ), and the 
contribution of the "new" dangerous terms vanish. This finishes the proof. 

E. Double soft limit of Goldstone boson amplitudes 

In this appendix we will discuss the properties of the on-shell scattering amplitudes of the Goldstone 
bosons, which are dictated by the symmetry, namely the limits of the amplitudes for soft external 
momenta. Some of these properties have been obtained in the special case of pions by PCAC 
methods in the late sixties (see e.g. [48]). Here we enlarge and reformulate them in a more general 
form appropriate for our purposes with stress on the proof of the double soft limit discussed recently 
for pions and A^ = 8 supergravity in [50]. 
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Let us assume a general theory with spontaneous symmetry breaking according to the pattern 
G —)■ H where the homogeneous space G/H is a symmetric space, i.e. the vacuum little group H 
is the maximal subgroup invariant with respect to some involutive automorphism of G ( "parity" ) . 
This implies the following structure of the Lie algebra of G 



(E.l) 



[T a ,T h ] 


= ifft c T c 


[T a ,X b ] 


= ifx bc X c 


[X a ,X b ] 


■ rpabcn-ic 



Here T a and X a are the unbroken and broken generators respectively and /<p , fx and F are 

fabt 
'X 



the structure constants. The chiral nonlinear sigma model is a special case for which /j. = f% 



rpabc fabc 

The invariance of the theory with respect to the group G can be expressed in terms of the 
Ward identities for the correlators in the general form 

n 

ip fl (K(p)d 1 ( Pl )...d n ( Pn )) = -£i(6i(pi)...^6i(pi + p)...6„(p»)> (e.2) 

i=l 

n 

i P »(Ai(p)d 1 ( Pl )...d n ( Pn )) = -Y J ^Oi(pi)...5 x d l ( Pl +p)...d n ( Pn )}. (e.3) 

Here V®(x) and A®(x) are the Noether currents corresponding to the generators T a and X a re- 
spectively (in analogy with the chiral theories we will call them vector and axial currents in what 
follows and to the Ward identities (E.2) and (E.3) we will refer to the vector and axial WI) , Oi(x) 
are (generally composite) local operators, S^Oi(x) and 5 x Oi(x) are their infinitesimal transforms 
with respect to the generators T a and X a . The tilde means the Fourier transform 

Oi(p)= [ d 4 xe [ P- x Oi(x). (E.4) 

According to the Goldstone theorem the spectrum of the theory contains as many Goldstone 
bosons ir a as the broken generators X a for which the currents A^x) play the role of the interpolating 
fields, i.e. 

(0\A;(0)\ir b (p)) = ip,F5 ab . (E.5) 

where F is the Goldstone boson decay constant. Let as denote M ai '" a "(pi, . . . ,p n ) the on-shell 
scattering amplitude of the Goldstone bosons ir ai (pi), . . . ,ir an (p n ). In what follows we will con- 
centrate on the properties of M ai,,,an (pi, . . . ,p n ) dictated by the symmetry, i.e. those which are 
encoded in the WI (E.2) and (E.3). 

E.l Vector WI and symmetry with respect to H 

The invariance with respect to the unbroken subgroup H implies 



n 

\ "* j;Cia 

i=l 



jaa i o M a 1 ...a i . 1 ba i+1 ...a n(<pu . . . jPn ) = Q. (E.6) 
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This can be understood as the consequence of the vector WI of the form 



-ipf* 



<V£(pM£(pi) • • • Ali( Pn )) = -J2i(A°\(pi) ■ ■ ■ 5 a T A%{p + Pi) ■ ■ ■ A^iPn)) (E.7) 



j=l 



Note that the infinitesimal transformations 5 a V b and S a A^ of these currents with respect to the 
generator T a of the unbroken subgroup H are as follows 

%Al = -ftK (E.8) 

W = -fr bc V u c . (E.9) 

Because there is no pole for p — > in the correlator on the left hand side of (E.7), we get in this 
limit 



E/ aai6 W 1 (Pi)---^(R)-^(Pr.)> =0. (E.10) 

i=l 

Using the LSZ formula we get according to (E.5) 

(A;\( Pl )...A%(p n )) = [\{\zA M^- a -( Pl ,...,p n ) + R^;- (E.ll) 

where Z^ = iFp^ and the remnant R^ " is regular on shell in the sense that 

iim (n^)^i::: = o. (E.12) 



pf^O 



\i=l 



which implies (E.6) for the on-shell amplitude M 01- " a ™(pi, . . . ,p n ). 

E.2 Soft vector current singularity 

Let us assume now the following matrix element 

<V£(p)|7T ai (Pi) . . . T^fa) • • • 7T a "(p„)>. (E.13) 

In what follows we will discuss the behavior of this object in the limit p — > 0. On the level of the 
Feynman graphs, the only singularities in the soft limit p — > are those which stem from the one- 
Goldstone-boson-reducible graphs for which the vector current V£(p) is attached to the external 
Goldstone boson line. The potential singularities are therefore of the form (see Fig. 14 ) 

(V^ P )^(0)\^(p 1 )) 1PI iA a ^((p - Pi ) 2 )(r k (0)K 1 (Pl) ■ ■ ■ T^iPi) ■ ■ ■ K an {Pn))lPI (E.14) 

where the subscript 1PI means one-Goldstone-boson-irreducible block, the hat means omitting of 
the corresponding particle, (j) a {x) is the Goldstone boson interpolating field normalized as 

(O|0 a (O)|7r 6 (p)) = 5 ab (E.15) 
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Pi, ai Pi-!, <Zj_i 

P, a, fi 

/ 




^_^_^_ Pi, a t 



v n . a n Dij.1, a,- j-i 



Figure 14: Graphical representation of the singular contributions to the matrix element (E.13). 

and A ajak (q 2 ) is a Goldstone boson propagator. For q 2 — > we have 

A a ^(q 2 ) = ^(l + 0(q 2 )). (E.16) 

As a consequence of the Lorentz invariance, invariance with respect to H and LSZ formulae we 
have 

(V*(p)^(0)K>( Pl )) 1PI = if™ iaj F v (p 2 )(2 Pi - p)„ + 0((p - Pi ) 2 ) (E.17) 

where Fy(p 2 ) is the on-shell vector form-factor defined as 18 

(^(p-Pi)\V-(p)\^( Pi )) = if a x a ^F v (p 2 )(2 Pl -p),. (E.18) 

We can fix the normalization of the vector currents V^ in such a way that 

F v (p 2 ) = l + 0(p 2 ). (E.19) 

Analogously we have 

(0 afc (O)|7T ai (Pl) • • • T^Pi) ■ ■ ■ K a "(p n ))lPI = M ai '" a '- 10fca *+ 1 - a »(p 1 , . . . , Pn ) + 0{(p - Pi f). (E.20) 

Using {p — Pi) = —2{p ■ p^ + p 2 and putting all the ingredients together we get for p — > 
<^»K ai (pi) . . . vt^(k) . . . Tr a "(p n )) = J2 ^^M«"^^^(pi, . . . ,p n ) + O(l) 

■ =1 ^J>-Pij 

(E.21) 

E.3 Axial WI and Adler zero 

To illustrate the method which we will use in the next subsection, let us briefly recapitulate the 
textbook example of the derivation of the Adler zero for the amplitude M ai '" an (pi, . . . ,p n ) (see 
e.g. [51]). Let us start with the axial WI in the form 

n 

-i^(2>)l»i(p 1 )...l^(p n )) = -^i{^(pi)--^i^(p + Pi)---^(Pn)) (E.22) 



3 The form of the right hand side is dictated by H-invariance, Bose and crossing symmetry. 
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where now 

5 a x A h v = -F abc VZ 

SxV„ b = -ff c K- (E.23) 

Applying on both sides of (E.22) the LSZ reduction to all but one axial currents, we get the 
conservation of the axial current in terms of the transversality of the matrix element of A between 
the initial and final states \i) and (/| 

-ip^/|2»|*>=0. (E.24) 

On the other hand from (E.ll) we get the Goldstone boson pole dominance for p 2 — > 

- W|1»K) = XpTZ^f + n a (p)\i) ~ ^K,fi ( E - 25 ) 

where Z^ = iFp^ and the remnant R a fi is regular in this limit 

lim p 2 R a n fi = 0. (E.26) 

Putting (E.24) and (E.25) together we get for the amplitude with emition of the Goldstone boson 
ir a (p) in the final state 

{f + 7r a (p)\ l ) = ^R^ fl . (E.27) 

Provided the following stronger regularity condition holds 

limp^/i = 0, (E.28) 

we get 

</ + ^(0)|i>=0, (E.29) 

i.e. the Adler zero for p — > 0. 

An useful off-shell generalization of the formula (E.25) reads 

-ip"(A£(p)l£ (px) . . . AH(jp n )) = iF(n a (p)\A;\ ( Pl ) . . . AH(p n )) - ^R a ^;:. (E.30) 

where 

lim p 2 R a '2- = 0. (E.31) 

p2^ M ' M1 - 

and using the Ward identity (E.22) and (E.23) we get 
F{n a (p)\A*\( Pl )...A° : (p n )) 

n 

= P ^RZ\::. + Y, FaaiC Wp^---K(p+p^--- I Z(pn))- (e.32) 



52 



E.4 Double soft limit 

Our starting point is the axial WI (E.22) rewritten in the form 

-^(A^p)Al( q )A^( Pl )...A^(p n )) = -i{8\Al(jp + q )Al\(p{)...All(jp n )) (E.33) 

n 

i=i 
Multiplying then both sides by — \q v and using the axial WI (E.22) once again we get 

= q u F abc (V u c (p + q)Al\{ Pl ) . . . A*l(p n )) 

n 

+ £ F^F^ d {Al\{p 1 )...V^{p l + q )...V^{p + p 3 )...All{ Pn )) 

n 

+ £ F a ^fY{Al\ (px) . . . A% (p + q + th) . . . AH ( Pn )). (E.34) 

The left hand side of (E.34) is symmetric with respect to the interchange of (p, a) -H- (q, b); its right 
hand side can be therefore rewritten in the manifestly symmetric form 

-p^q u {Al{p)Al{q)Al\{ Pl ) . . . A-l(p n )) 
= ~\{P ~ qYF ahc (V u c (p + q)A«\( Pl ) . . . A%( Pn )) 

n 

+ E F aa ^F^ d (Al\ ( Pl )... V£ (pi + q)... V£. (p + Pj )... A% (p n )) 

+ \JZ (F aa * c f b x cd + F ba > c ff d ) (A% (px) ...A%(p + q + Pi )... A%(p n )). (E.35) 

i=i 

On the other hand, the LSZ formula gives for p 2 , q 2 — > 

(Al(p)Ai(q)Al\( Pl )...A^( Pn )) = Y J ^(0\Al\n c ( P ))^(0\A b / \n d (q)) 

c.d 



x ^{p)**{q)\A% ( Pl ) . . . Afcfa)) + R$- (E.36) 
where the regular remnant satisfies 



lim pV<,'-=0. (E.37) 



Therefore, using (E.5) we get 



-pV<^>)^Ki(Pi) : ..^:(Pn» 

F 2 {n a (p)Aq)\Al\ (Pi) • • • A«i( Pn )) - p»q"Rf>- (E.35 
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On the other hand applying the LSZ reduction to (E.34, E.35) (let us note that only the first terms 
on the right hand side has the appropriate poles at p 2 ,q 2 — > 0) we get 

P»q u (A*(p)A b u ( q )\7r a Hpi) • • • APi) ■ ■ ■ ^ n (Pn)) 
= q»F abc (K(p + q)\^ ( Pl ) . . . APi) ■ ■ ■ v an {Pn)) 
= p^F bac (V^(p + q)\^ ( Pl ) . . . APi) ■■■* a " (Pn)) 

= -\F ahc {p - qyiyfa + q)\^ (pi) . . . n d (p,) . . . vr a "(p n )> (E.39) 

and as a consequence of LSZ reduction of (E.38) 

F 2 (n a (p)7v\q)\^( Pl ) . . . 7r ai ( Pi ) . . . vr a "(p n )> 
= -\F ahc {p - qT(Vfe + q)\7T a H Pl ) . . . 7T d ( Pi ) . . . K an (p n ))+P»q u R$-\LSZ. (E.40) 

According to (E.21) we have for p, q — >■ 

-\F ahc {p - qTiVfa + q)\^{pi) . . . n d &) ■ ■ ■ ^(Pn)) 
= ~ E FabC fT d {2Pl ~^ q) qy ^^\ ^ (Pl)--- Ap>) • • — ( Pn )) + 0(p - q) 



1 J- F ^ fT d P_i • (P g) ^ {pi) nd{pi) n a n {pn)) 



2^ X Pi-(p + q) 

+0 (p - q, **'?. ) (E.41) 

V Pi- (p + q)J 

For p 2 = q 2 = we finally get 

FoV^VWK 1 (Pl) ■ ■ ■ * ai (Pi) ■ ■ ■ * an (Pn)) 
= ~\JZ F^fT ^ '^^l (* ai (Pl) ■ ■ ■ APi) ■ ■ ■ *"» (Pn)) 

2 ~[ Pi- (P + Q) 

+p^R$-\ LSZ+ 0(p-q). (E.42) 

Provided condition stronger than (E.37) holds, namely lim p :q ^o p^q u R^ '"\lsz = (cf. (E.28)), we 
get as a result 

limF 2 (7r a (tp)7r 6 (ig)|7r ai (pi) . . . vr^(p 4 ) . . . vr a "(p n )> 

= ~\ E F^ ^f"^ /' (Pi) • • • APi) . . . vr a "(p n )>. (E.43) 

2 ~y Pi' (p+q) 

For the chiral nonlinear sigma model corresponding to the symmetry breaking G x G — > G, we have 
pabc _ jahc _ jabc an( j we g e ^. ^ e formula (5.3) as a special case. 
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